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1.0  SUMMARY 


A  supersonic  nozzle  with  supersonic  iodine  injection  was  designed  and  studied  with 
Planar  Laser-Induced  Fluorescence  (PLIF).  The  nozzle  simulates  Chemical  Oxygen  Iodine 
Laser  (COIL)  flow  conditions  with  non-reacting,  cold  flows.  A  laser  sheet  near  565  mn  excited 
the  iodine,  and  the  fluorescence  was  imaged  with  an  intensified  and  gated  CCD  camera. 
Streamwise  and  semi-spanwise  (oblique -view)  images  were  taken,  where  the  presence  of 
injected  flow  was  highlighted.  With  these  images,  the  flow  structures  were  identifiable  and  the 
mixing  quality  between  the  primary  and  injected  flow  was  quantitatively  measured  with 
histograms,  structure  size  measurements,  and  jet  penetration. 

Four  different  injection  scenarios  are  investigated.  The  first  scenario  includes  a  single 
injector  positioned  downstream  of  the  nozzle  throat.  To  enhance  the  mixing  between  the  flows, 
trip  jets  were  placed  in  the  wake  of  the  single  jet.  The  trip  jets,  significantly  smaller  than  the 
primary  iodine  jet,  are  intended  to  destabilize  the  counter-rotating  vortex  pair  (CRVP)  of  the 
primary  jet.  Three  different  trip  jet  configurations  are  compared  for  their  ability  to  enhance 
mixing  between  the  oxygen  and  iodine  flows.  While  the  mixedness  between  the  injected  and 
primary  flows  was  not  optimized,  the  comparison  between  the  four  different  injection  scenarios 
may  aid  in  future  nozzle  design  and  mixing  optimization. 
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2.0  INTRODUCTION 


2.1  Chemical  Oxygen  Iodine  Lasers 

In  1977,  the  Air  Force  Weapons  Laboratory  demonstrated  the  Chemical  Oxygen  Iodine 
Laser  (COIL).1  McDermott  and  coworkers  passed  chlorine  gas  vapor  through  a  rotating  disk 
device,  which  contained  basic  hydrogen  peroxide  as  a  fuel,  to  obtain  the  excited  state  of 
molecular  oxygen,  O2LA).  Molecular  iodine  was  injected  into  the  gaseous  flow  of  excited 
oxygen,  02(  A)  and  02(  2),  and  ground-state  oxygen,  02(  2),  to  produce  the  excited  state  of 
atomic  iodine.  The  “P1/2  -  P3/2  transition  of  atomic  iodine  produces  an  output  wavelength  of 
1.315  pm.  The  whole  transition  is  pumped  by  the  singlet  delta  excited  state  of  molecular  oxygen 
CL^A),  triggering  the  following  sequence  of  reaction  equations: 

02('a)+I2  — >02(32)+/2 

0!(,e)+/'^0!(!e)+2/(iP:!,2)  (1) 

Truesdell  et  al.,  Avizonis  and  Neumann,  and  Avizonis  and  Truesdell  give  detailed  descriptions  of 
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the  COIL  chemistry  and  mechanics.  ’  ‘ 

The  gases  are  typically  mixed  in  a  supersonic  nozzle,  with  excited  oxygen  and  a  diluent 
(usually  helium)  in  the  primary  flow  and  iodine  and  a  diluent  (helium)  in  the  injected  flow. 
Supersonic  nozzles  are  used  to  maximize  the  temperature  drop  downstream  of  the  throat  of  the 
nozzle  (and  thus  move  the  state  of  the  flow  farther  away  from  thermodynamic  equilibrium, 
facilitating  lasing)  and  to  increase  the  flow  per  unit  nozzle  throat  area.  The  low  temperature  of  a 
supersonic  nozzle  optimizes  the  production  of  I  (  P1/2).  Historically,  the  iodine  and  its  diluent 
have  been  injected  into  the  subsonic  region  of  the  nozzle  upstream  of  the  throat.  Injecting  the 
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iodine  in  the  subsonic  region  promotes  mixing  between  the  injected  iodine  and  the  primary 
oxygen  flows,  whereas  diffusive  mixing  is  difficult  in  the  supersonic  region  because  of 
compressibility  effects.  A  well-mixed  medium  increases  the  interaction  between  02(  'A)  and 
iodine,  to  produce  more  I  (2Pi/2).  A  well-mixed  medium  also  produces  a  more  uniform  laser 
beam  profile.  Numerical  studies  have  been  performed  for  chemical  lasers  with  iodine  injection 
in  the  subsonic  region  of  the  nozzle  by  Madden  and  Miller,  Miller  et  ah,  Madden  and  Solomon, 
and  Masuda  et  al.5'9  These  simulations  were  perfonned  to  investigate  the  spatial  and  temporal 
gradients  in  the  flowfield,  which  impact  the  laser  beam  quality. 

While  the  injection  of  iodine  has  traditionally  been  in  the  subsonic  region  of  the  nozzle, 
there  are  advantages  to  injecting  the  iodine  in  the  supersonic  region  of  the  nozzle.  Designs  in 
which  iodine  is  injected  in  the  supersonic  region  have  the  potential  to  provide  improved 
performance.  Injecting  the  iodine  downstream  of  the  nozzle  throat  decouples  the  injected  flow 
from  the  Singlet  Oxygen  Generator  (SOG),  the  device  that  chemically  makes  and  delivers 
02(1A).  This  decoupling  allows  changes  to  be  made  to  the  injected  mass  flow  without  affecting 
the  pressure  upstream  in  the  flow  from  the  SOG.  When  injected  into  the  subsonic  region,  the 
iodine  and  diluent  gases  are  a  significant  portion  of  the  SOG  molar  flow.  By  injecting  the  iodine 
and  diluent  helium  into  the  supersonic  region,  the  pressure  increase  and  velocity  decrease 
associated  with  injection  are  avoided,  reducing  the  loss  of  CL^A)  to  quenching  processes.  Lower 
pressure  and  increased  velocity  decrease  the  rate  of  CL^A)  loss,  and  a  greater  velocity  decreases 
the  time  it  takes  to  transport  CL^A)  downstream  of  the  throat,  providing  more  Chi'A) 
downstream  for  laser  power  extraction. 

Injecting  the  iodine  in  the  supersonic  region  does  hinder  the  mixing  between  the  primary 
flow  of  oxygen  and  the  injected  iodine  flow.  Mixing  between  supersonic  streams  is  more 
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difficult  because  of  the  increased  flow  stability  associated  with  compressible  flows  and  shock 
losses  that  decrease  jet  penetration.  The  supersonic  nozzle  and  the  iodine  injectors  must  be 
designed  to  provide  adequate  penetration  of  the  iodine  stream  into  the  oxygen  flow  and  adequate 
mixing  between  the  oxygen  and  iodine  flows  to  completely  dissociate  the  iodine  molecules.10 
The  process  of  dissociation  of  the  iodine  molecules  by  CL^A)  is  sensitive  to  the  rate  of  mixing. 
If  the  local  molecular  iodine  concentration  drops  too  rapidly,  dissociation  is  decreased  and  laser 
performance  drops.  COIL  nozzle  designs  with  supersonic  injection  of  iodine  have  been  studied 
numerically  by  Madden  and  coworkers  and  experimentally  by  Rosenwaks  and  coworkers  and 
Barmashenko  and  coworkers.10'12 


2.2  Fluid  Mechanics  Related  to  COIL  Nozzles 

The  flow  physics  of  iodine  injection  into  the  oxygen  flow  of  a  COIL  nozzle  are  very 
similar  to  those  of  a  fundamental  problem  of  a  jet  in  crossflow  (JICF).  There  have  been 
extensive  studies  of  JICF  structures  and  methods  used  to  perturb  the  stability  of  the  dominant 
structures  to  enhance  mixing.  This  research  and  methods  of  perturbation  are  directly  applicable 
to  COIL  nozzles,  with  a  caveat  that  COIL  nozzles  operate  under  relatively  low  pressures  when 
compared  to  many  other  flow  regimes. 

2.2.1  Jet  in  Crossflow  Structure 

Fric  and  Roshko  give  a  thorough  description  of  four  structures  present  when  an 

13 

incompressible  jet  enters  a  crossflow,  with  emphasis  on  the  structure  in  the  wake  of  the  jet. 

They  experimentally  studied  the  jet  in  crossflow  with  smoke- wire  visualization,  which  places 

closely-spaced  streaklines  into  the  flow.  In  their  experiments,  the  crossflow  Reynolds 
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number  ranged  from  3800  to  1 1,400  and  the  jet-to-crossflow  velocity  ratio  ranged  from  2  to 
10.  They  refuted  the  common  assumption  that  there  is  a  strong  similarity  between  the  near 
wake  of  a  jet  and  that  of  a  circular  cylinder.  In  their  experiments,  they  identified  four  types 
of  coherent  structures  in  the  near  field  of  the  jet:  1)  jet  shear-layer  vortices;  2)  system  of 
horseshoe  vortices;  3)  counter-rotating  vortex  pair  (CRVP);  and  4)  wake  vortices.  These 
structures  are  illustrated  in  Figure  1.  Initially,  the  jet  shear-layer  vortices  are  the  dominant 
flow  structures.  They  result  from  the  Kelvin-Helmholtz  instability  of  the  annular  shear  layer 
that  lies  on  the  leading  edge  of  the  jet  orifice.  The  horseshoe  vortices  arise  from  the  adverse 
pressure  gradient  upstream  of  the  jet  interacting  with  the  approaching  low-momentum  fluid 
in  the  wall  boundary  layer.  The  CRVP  is  the  dominant  system  of  primary-secondary  flow 
interaction  and  persists  far  downstream.  In  the  development  of  the  CRVP,  the  crossflow 
skirts  laterally  around  the  jet,  shearing  the  jet  fluid  along  its  edges  and  then  folding  the  face 
of  the  jet  over  itself  to  form  the  CRVP.  Finally,  the  wake  vortices  are  formed  in  the  region 
downstream  of  the  jet  (between  the  jet  and  the  wall).  Each  wake  vortex  has  a  termination  on 
the  wall  and  on  the  jet,  whereas  with  a  circular  cylinder  the  wake  vortices  terminate  only  on 
the  wall. 

Fric  and  Roshko  noted  that  it  was  essential  to  establish  the  source  of  vorticity  of  the 
wake  structures  to  understand  them.  From  their  smoke  visualization,  they  found  that  the 
vortical  structures  in  the  wake  of  a  jet  are  formed  from  the  vorticity  originating  in  the 
boundary  layer  of  the  crossflow  wall.  Alternating  separation  events  occur  on  each  side  of  the 
jet,  on  the  crossflow  wall,  resulting  in  eruptions  of  boundary  layer  fluid  and  vorticity  that 
interact  with  the  jet  fluid.  Fric  and  Roshko  conclude  that  the  wake  vortices  may  affect 
unsteady  components  of  the  CRVP  structure. 


5 


(''•unii-r-rnu'iny 

full 


Figure  1.  Cartoon  depicting  the  four  types  of  vortical  structures  associated  with  a  jet-in-crossflow.13 

Rivero  et  al.  also  performed  an  experimental  study  to  identify  the  incompressible 
JICF  structures.14  They  performed  hot-wire  measurements  and  PLIF  flow  visualization  on  a 
JICF  with  a  velocity  ratio  of  3.8  and  crossflow  Reynolds  number  of  6600.  They  concluded 
that  three  structures  (folded  vortex  rings,  horseshoe  vortices,  and  handle-like  structures) 
contribute  to  a  coherent  velocity  field.  The  folded  vortex  rings  result  from  the  roll-up 
process  of  the  shear  layer,  as  in  the  shear-layer  vortices  of  Fric  and  Roshko’s  work.  The 
handle-like  structures  link  the  boundary  layer  vorticity  with  the  CRVP  through  the  upright 
wake  vortices.  The  CRVP  is  strongly  modulated  by  this  coherent  velocity  field.  Their 
results  indicate  that  the  CRVP  is  not  a  steady  feature  of  the  jet  in  crossflow,  but  the  intensity 
of  each  counter-rotating  vortex  has  fluctuations  driven  by  Kannan-like  vortices  present  in  the 
jet  wake. 

Though  the  studies  by  Fric  and  Roshko  and  Rivero  et  al.  detail  the  structures  of  a  jet 
issuing  into  a  subsonic  crossflow,  the  description  of  the  jet  structures  still  applies  to 
supersonic  flows.13’14  Gruber  et  al.,  in  their  studies  of  the  compressibility  effects  of  a 
supersonic  JICF,  developed  a  three-dimensional  schematic  of  the  flow  structures,  as  shown  in 
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Figure  2. 15  Because  the  crossflow  is  supersonic,  a  bow  shock  forms  just  upstream  of  the  jet, 
increasing  the  jet  stability.  The  jet  is  supersonic  and  forms  a  shock  structure  (Mach  disk) 
downstream  of  the  jet  exit.  The  wake  region,  horseshoe  vortex  region,  and  the  CRVP  are  all 
present  with  a  supersonic  jet  issuing  into  a  supersonic  crossflow. 


Figure  2.  A  three-dimensional  schematic  of  the  structures  in  a  compressible  (supersonic)  JICF 
system.15 

2.2.2  Perturbation  of  Injected  Flow 

The  jet  structures,  mainly  the  CRVP,  can  be  perturbed  to  produce  unsteadiness  or 

altered  in  shape  to  increase  the  mixing  between  the  jet  and  crossflow.  Voignier  et  al.  studied 

three  different  methods  to  perturb  the  jet  structure,  motivated  by  the  need  to  stretch  the 

mixing  interface  in  a  deuterium  fluoride  (DF)  combustion  driven  chemical  laser  and  increase 

mixing.16  They  compared  small  signal  gain  profiles  from  different  perturbation  methods: 

trip  jet  nozzles,  ramp  nozzles,  and  deflector  nozzles.  The  trip  jets  used  either  inert  gas  or 

cavity  fuel  and  were  injected  transversely  between  reactant  streams  through  small  orifices. 

In  that  study,  the  metric  for  comparing  the  perfonnance  of  the  trip  jets  to  the  ramp  and 

deflector  nozzles  involved  both  perfonnance  and  manufacturing  costs.  To  compare  the 
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mixing  enhancement  capability,  the  small  signal  gain  and  the  distance  of  positive  gain  after 
the  nozzle  exit  plane  (NEP)  were  measured.  With  no  trip  jets,  ramps,  or  deflectors,  the 
maximum  small  signal  gain  was  1.5%/cm  and  the  distance  of  positive  gain  was  40  mm  from 
the  NEP.  The  nozzle  with  the  trip  jets  produced  the  largest  small  signal  gain,  with  a 
maximum  of  7%/cm,  with  positive  gain  measured  for  over  30  mm  from  the  NEP.  The 
deflector  and  ramp  nozzles  had  lower  small  signal  gain  and  the  region  of  positive  gain  was 
less  than  20  mm.  With  the  trip  jet  nozzle  having  a  higher  peak  gain  and  longer  region  of 
positive  gain,  it  was  concluded  that  the  mixing  rate  was  higher  than  that  of  the  ramp  and 
deflector  nozzles. 

Cenkner  and  Driscoll  used  Planar  Laser-Induced  Fluorescence  (PLIF)  and  Laser 
Doppler  Velocimetry  (LDV)  to  study  supersonic  cold-flow  mixing  of  nozzles  that  employ 
gas  trips.17  Their  nozzle  simulated  chemical  laser  flow  conditions,  with  the  reacting  species 
being  injected  adjacently  in  a  nozzle.  The  trip  jets  were  positioned  near  the  NEP  of  both 
nozzles  of  the  reacting  species  and  utilized  nonreacting  flow.  Cenkner  and  Driscoll  conclude 
that  there  are  three  convective  mechanisms  that  could  contribute  to  the  enhanced  mixing  of 
adjacent  gas  streams:  1)  a  low-pressure  wake  of  the  trip  jets  that  draws  gas  from  the  adjacent 
flow;  2)  a  high-pressure  region  that  is  established  in  the  adjacent  flow,  that  pushes  this  flow 
into  the  main  flow;  and  3)  a  trip-jet  induced  breakup  that  possibly  triggers  a  spreading  effect 
in  the  adjacent  flow.  They  also  concluded  that  trip  jets  did  not  increase  the  maximum 
turbulence  intensity,  but  instead  spread  the  low-level  turbulence  in  the  wake  of  the  base 
region  into  the  core  of  the  nozzle  flows. 

Driscoll  published  experimental  results  comparing  ramp  and  trip  jet  nozzles 
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applicable  to  DF  chemical  laser  nozzles.  The  trip  jets  supplied  inert  gas  into  an  inert  gas 
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layer  separating  the  reactant  streams,  retarding  the  molecular  mixing.  The  inert  gas  also 
perturbed  the  density  in  the  laser  cavity,  potentially  degrading  the  laser  beam  uniformity. 
Driscoll’s  LIF  (laser-induced  fluorescence)  experimental  results  indicate  that  the  ramp  and 
trip  jet  nozzles  (with  inert  gas)  produce  similar  levels  of  mixing  enhancement.  The  ramp  and 
trip  nozzles  both  doubled  the  reactant  interface  length  within  a  centimeter  of  the  NEP. 
However,  because  of  the  possible  inert  gas  layer  fonned  by  the  trip  jets  in  that  case,  in  future 
work  it  may  be  advantageous  to  use  a  reactant  gas  in  the  trip  jet. 

While  Voigneir  et  al.,  Cenkner  and  Driscoll,  and  Driscoll  employed  trip  jets  to 
destabilize  both  the  primary  and  secondary  flows,  Takeuchi  et  al.  employed  trip  jets  with  a 
reactant  gas  to  destabilize  the  secondary  flow.16'19  An  ejector  COIL  consisting  of  two- 
dimensional  slit  nozzles  with  trip  jet  iodine  injection  was  designed  and  tested.  The  trip  jets 
were  located  at  the  NEP  of  the  nitrogen  nozzle  (which  was  adjacent  to  the  oxygen  nozzle). 
The  iodine  was  mixed  into  the  nitrogen  flow  via  the  trip  jets  and  then  the  nitrogen/iodine 
flow  was  mixed  into  the  adjacent  oxygen  flow.  They  used  a  horizontal  Pitot  scan  to 
experimentally  demonstrate  that  the  mixing  ability  of  the  trip  jets  was  sufficient.  While  the 
slit  nozzle  design  is  not  similar  to  the  supersonic  injection  scheme  studied  here,  it  does 
demonstrate  the  mixing  enhancement  of  trip  jets,  because  they  not  only  injected  the  iodine, 
but  also  were  able  to  enhance  the  mixing  between  the  nitrogen  and  oxygen  flows. 

The  injected  flow  may  also  be  altered  by  means  other  than  ramps,  deflectors,  and  trip 
jets.  The  jet  may  be  influenced  by  forcing  the  injected  flow  at  a  selected  frequency  or 
imparting  the  jet  with  a  swirling  component  before  it  issues  into  the  cross  stream  flow.  A 
forcing  frequency  acts  upon  existing  instabilities  in  the  flow,  whereas  a  swirling  jet 
introduces  a  new  instability.  Murugappan  et  al.  used  Rayleigh/Mie  scattering  from  naturally 
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occurring  ice  crystals  to  study  the  mixing  and  penetration  of  a  forced  jet  in  a  Mach  2  cross 
stream.  They  forced  the  underexpanded  jet  into  the  cross  stream  with  a  high-frequency 
actuator  at  three  different  frequencies  (0  Hz,  900  Hz,  and  5  kHz).  They  found  that  both  of 
the  forced  cases  provided  superior  mixing  characteristics  and  higher  penetration  of  large- 
scale  structures  into  the  freestream  than  the  unforced  case.  In  a  subsequent  study, 
Murugappan  et  al.  used  Rayleigh/Mie  scattering  and  NO  (nitric  oxide)  PLIF  to  study  the 
mixing  and  penetration  characteristics  of  a  swirling  injector  that  consisted  of  two 
independent  streams.”  One  stream  was  annular  with  swirl  vanes  at  the  exit  to  impart  swirl 
on  the  jet.  The  other  stream  was  a  central  jet  to  control  the  penetration  of  the  jet  into  the 
freestream  and  prevent  early  vortex  breakdown.  The  optimal  case  that  was  tested  showed  an 
increase  of  16%  in  mixing  area  with  the  streamwise  images  and  an  increase  of  78%  in  total 
area  contained  within  the  jet  boundary  with  the  cross-stream  images. 

The  shape  of  the  injector  also  affects  penetration  and  large-scale  structure  size  of  the 
jet  in  crossflow.  Gruber  et  al.  and  Gruber  et  al.  used  planar  Rayleigh/Mie  scattering  to  study 
two  sonic  transverse  injectors.  '  The  experiment  of  Gruber  et  al.,  involving  a  circular 
injector,  was  simulated,  with  comparable  results,  using  the  Reynolds  Averaged  Navier  Stokes 
(RANS)  equations,  by  Palekar  et  al."  ’  Of  the  two  sonic  transverse  injectors  studied  by 
Gruber  et  al.,  one  injector  was  circular  and  the  other  was  elliptical  (with  the  major  axis 
aligned  in  the  lateral  direction)."  Both  injectors  were  directed  perpendicular  to  the 
crossflow.  The  elliptical  jet  was  found  to  spread  more  rapidly  in  the  lateral  direction  than  the 
circular  jet,  indicating  an  axis  switching  phenomenon  within  the  jet.  However,  the  elliptical 
jet  did  suffer  a  20%  reduction  in  penetration  compared  to  the  circular  jet.  The  elliptical  jet 
was  wide  and  flat  with  the  counter-rotating  vortices  spread  farther  apart  than  the  vortices 
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from  the  circular  jet.  This  spread  in  vortices  provided  a  larger  central  entrainment  region 
where  freestream  fluid  may  be  entrained  into  the  jet  core.  Though  the  penetration  of  the  jet 
suffered,  the  increased  entrainment  region  enhanced  the  mixing  of  the  jet  and  crossflow. 
Blanchard  et  al.  and  Landman  and  Saffman  discuss  the  nature  of  a  JICF  using  elliptical 
geometry.25'26  They  conclude  that  the  elliptical  cross  section  of  a  CRVP,  regardless  of  orifice 
shape,  can  cause  the  loss  of  jet  stability.  This  loss  in  stability  would  also  have  the  potential 
to  enhance  the  mixing  between  the  jet  and  crossflow. 

The  work  by  Viti  et  al.  examines  the  use  of  smaller  jets  to  alter  the  flow  behind  a 
primary  jet,  though  it  does  not  elaborate  on  the  structural  interactions  between  a  primary  jet 
and  an  array  of  smaller  jets.  Viti  et  al.  numerically  simulated  and  experimentally  studied  a 
primary  jet  issuing  into  a  Mach  4  crossflow  with  an  array  of  smaller  injectors  downstream  of 
the  primary.  The  jets  were  all  sonic  and  injected  perpendicular  to  the  crossflow,  with  the 
smaller  injectors  injecting  0.5  to  1%  of  the  mass  flow  rate  of  the  primary  jet.  They  used 
pressure  sensitive  paint  to  measure  the  pressure  field  and  Schlieren  photographs  for  imaging. 
The  numerical  simulation  was  done  with  RANS  equations  and  the  Wilcox  k- co  turbulence 
model.  They  did  not  measure  the  impact  the  smaller  jets  had  on  fluid  mixing,  but  did  find 
that  the  smaller  injectors  reduced  the  size  and  intensity  of  the  low-pressure  region 
downstream  of  the  primary  jet. 

2.3  COIL  Diagnostics 

Nonintrusive  imaging  of  the  interaction  between  the  oxygen  and  iodine  flows  can  be  used 
to  determine  the  degree  of  mixing  between  the  primary  and  secondary  flows,  the  concentration 
of  the  iodine  flow,  and  the  temperature,  pressure,  and  velocity  of  the  flows.  Knowing  the 
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thennodynamic  characteristics  and  the  mixing  quality  between  the  oxygen  and  iodine  flows  can 
be  useful  in  optimizing  the  COIL  mixing  nozzle. 

There  are  several  different  types  of  flow  imaging,  including  Particle  Image  Velocimetry 
(PIV),  Rayleigh/Mie  scattering,  Schlieren  imaging,  and  Planar  Laser-Induced  Fluorescence 
(PLIF).  As  an  example  of  PIV  in  a  similar  application,  Borg  et  al.  used  PIV  to  measure  particle 
velocity  in  a  low-pressure  turbulent  jet.”  The  jet  fluid  was  seeded  with  tracer  particles  that  were 
illuminated  by  a  plane  sheet  of  light.  A  charge-coupled  device  (CCD)  camera  acquired  images 
pairs  with  a  known  time  separation.  This  allowed  for  particle  displacement  mapping  from  one 
image  to  the  next.  To  make  such  mapping  possible,  the  time  separation  between  images  must  be 
small  enough  that  the  particle  movement  is  a  small  fraction  of  the  image  size. 

Rayleigh/Mie  scattering  is  similar  to  PIV  in  that  it  uses  a  laser  sheet  to  illuminate  the 
flow.  Rayleigh  scattering  is  observed  directly  off  tracer  gas  (e.g.,  biacetyl)  molecules,  whereas 
Mie  scattering  is  produced  by  light  scatter  from  tracer  particles  in  particle-seeded  flow.  A  CCD 
camera  is  used  to  capture  the  scattered  light.  The  image  intensities  are  directly  related  to  the 
tracer  concentration,  as  long  as  the  tracer  behaves  as  a  passive  scalar  during  the  mixing  process. 
An  example  is  the  work  of  Gruber  et  al.,  who  studied  JICF  flowfields  using  silicon  dioxide 
particles  as  tracers.15 

Schlieren  photography  is  also  used  to  image  density  gradients  in  the  flow  and  is  good  for 
imaging  the  boundary  of  a  JICF,  for  example.  Ben-Yakar  and  Hanson  used  fast-framing 
Schlieren  photography  to  study  the  time  evolution  of  JICF  structures.”  A  drawback  to  Schlieren 
imaging  is  that  while  it  can  be  used  to  identify  structures  along  the  edges  of  a  jet,  it  cannot 
discern  the  internal  features.  Another  limitation  is  that  it  integrates  the  density  gradients  along 
the  beam  propagation  path,  unlike  the  previously-discussed  methods  that  use  a  thin  sheet  of  light 
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to  illuminate  a  planar  section  of  the  flow,  and  integration  occurs  only  through  the  thickness  of 
the  light  sheet. 

PLIF  is  an  imaging  technique  where  a  thin  planar  sheet  of  laser  light  is  used  to  induce 
fluorescence  in  a  species  carried  by  the  flow.  The  fluorescent  light  (with  the  wavelength  distinct 
from  that  of  the  laser  light)  is  then  captured  with  a  CCD  camera.  The  laser  wavelength  used  is 
chosen  based  on  the  absorption  spectrum  of  the  fluorescing  species.  Some  common  fluorescing 
species  are  acetone  '  ,  nitric  oxide  ,  OH  for  reacting  flows  ,  and  iodine.  The  fluorescing 
species  is  chosen  based  on  ease  of  use,  level  of  toxicity,  or  whether  it  is  already  present  in  the 
system  to  be  imaged. 

Hiller  and  Hanson  discuss  the  use  of  iodine  in  PLIF  fluid  mechanic  experiments.34  They 
argue  that  saturating  the  fluorescence,  with  the  use  of  a  continuous-wave  excitation  laser,  allows 
the  fluorescence  signal  to  be  independent  of  the  quenching  rate  and  laser  intensity,  which  is 
desirable.  With  non-saturated  fluorescence  and  a  continuous-wave  excitation  laser,  the  PLIF 
signal,  Sf,  is  a  function  of  the  excitation  laser  power,  /,  fluorescence  decay  rate,  A2i,  collisional 
quenching  rate,  Q,  Boltzmann  population  fraction,  Fh  number  density  of  the  iodine  species  being 
irradiated,  Nn,  and  a  correction  for  the  optical  collection  geometry  and  molecular  transition,  C  ' 

Sf  =CI - - — F.N.  (2) 

This  direct  relationship  allows  the  determination  of  the  number  density  of  the  fluorescing  iodine 
molecules,  which  directly  relates  to  concentration  of  the  injectant. 

Iodine  PLIF  has  been  successfully  used  to  measure  temperature,  pressure,  and  velocity  in 
a  compressible  flow  field.  Hartfield  et  al.  used  the  fluorescence  broadband  to  determine 
temperature,  and  the  shape  and  position  of  the  fluorescence  spectrum  to  determine  velocity  and 
pressure  in  a  compressible  flowficld.  Their  uncertainties  were  quite  reasonable  with  6%  for 
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temperature,  5%  for  pressure,  and  25  m/s  for  velocity.  Iodine  PLIF  can  also  been  used  to 
measure  injectant  mole  fraction  and  iodine  concentration  in  a  flowfield.35,34 

With  iodine  present  in  COIL  systems,  using  PLIF  with  iodine  as  the  fluorescing  species 
is  an  obvious  choice  for  fluid  imaging.  There  are  several  benefits  to  using  iodine:  it  has  a  large 
absorption  cross  section,  it  fluoresces  strongly  in  the  visible,  it  has  a  high  vapor  pressure,  and  it 
has  numerous  transition  lines  across  a  large  range  of  frequencies  allowing  the  use  of  different 
laser  wavelengths  for  excitation. 


2.4  Thesis  Overview 

The  JICF  structure  analyses  of  Fric  and  Roshko  and  Rivero  et  al.  suggest  that  the  CRVP 
can  be  destabilized  or  stretched  by  producing  a  disturbance  in  the  wake  of  the  jet  that  mimics 
wake  vortices.13’14  The  wake  vortices  occur  in  the  downstream  region  of  the  jet  and  connect  the 
crossflow  boundary  layer  to  the  CRVP.  A  ramp,  deflector,  or  small  trip  jet  could  create  a 
disturbance  in  this  portion  of  the  flow  to  enhance  the  action  of  the  wake  vortices  or  create 
additional  wake  vortices  to  perturb  the  CRVP.  In  this  study,  trip  jet  interactions  with  a  single  jet 
will  be  investigated.  The  trip  jet(s)  will  be  located  in  the  subsonic  wake  region  of  the  jet  to 
connect  with  the  CRVP. 

The  trip  jet  designs  for  DF  and  COIL  chemical  lasers,  as  tested  by  Voignier  et  al., 
Cenkner  and  Driscoll,  Driscoll,  and  Takeuchi  et  al.,  demonstrated  the  capability  of  the  trip  jets  to 
enhance  mixing  between  two  components  of  the  flow.16’19  It  was  noted  that  using  reactants  in 
the  trip  jets  reduces  the  chance  of  having  an  inert  gas  layer  present  between  the  two  reacting 
species,  an  effect  that  would  inhibit  mixing  (and  thus  decrease  the  chemical  laser  efficiency). 
However,  trip  jets  were  not  the  only  tools  studied  to  enhance  mixing  with  injected  flow.  Though 
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trip  jets  seem  like  a  natural  device  to  stretch  the  structures  in  the  issuing  jets,  the  flow  may  be 
disturbed  or  influenced  by  forced  injection  of  the  jet  with  a  high-frequency  actuator,  or  swirl 
could  be  imparted  to  the  jet  at  its  exit.20'21 

The  objective  of  this  thesis  is  to  investigate  injection  penetration  and  mixedness  with 
supersonic  iodine  injection  in  the  supersonic  region  of  a  nozzle.  A  nozzle  was  designed  with 
COIL  chemistry  considerations  and  was  tested  using  non-reacting  flows.  A  single  jet  of  iodine 
and  diluent  helium  was  positioned  just  downstream  of  the  nozzle  throat  and  the  injected  iodine 
was  imaged  using  PLIF.  Using  the  PLIF  images,  the  jet  penetration  and  mixedness  of  the 
oxygen  and  iodine  flows  were  determined.  To  enhance  the  mixing  between  the  flows,  trip  jets 
were  placed  in  the  wake  of  the  single  jet.  The  trip  jets,  significantly  smaller  than  the  primary 
iodine  jet,  were  placed  with  the  intent  of  destabilizing  the  CRVP  of  the  primary  jet.  Different 
configurations  of  the  trip  jet(s)  are  compared  for  their  ability  to  enhance  mixing  between  the 
oxygen  and  iodine  flows. 

The  next  chapter  discusses  the  experimental  system  used  in  this  effort,  including  the  gas 
delivery  system  and  the  nozzle  design,  verification  of  the  nozzle  design,  and  the  specifications  of 
the  single  injector.  The  specifics  of  the  PLIF  diagnostics  are  detailed  in  Chapter  4.  Chapter  5 
presents  the  PLIF  results  with  the  single  injector,  while  Chapter  6  presents  the  results  of  the 
single  injector  with  three  different  trip  jet  scenarios.  Chapter  7  concludes  the  thesis  and  presents 
recommendations  for  future  work. 
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3.0  DESCRIPTION  OF  EXPERIMENTAL  SYSTEM 


3.1  Introduction 

This  chapter  describes  the  gas  delivery  system  and  the  nozzle  designed  for 
experimentation.  The  gas  delivery  system  includes  a  vacuum  pump,  delivery  of  the  primary 
gases  to  the  plenum  of  the  nozzle,  and  delivery  of  iodine  and  carrier  helium  to  the  injector(s)  of 
the  nozzle.  The  nozzle  was  designed  according  to  Chemical  Oxygen  Iodine  Laser  (COIL)  gas 
flow  specifications,  with  a  design  Mach  number  of  2.2.  The  nozzle  was  tested  to  verify  the 
Mach  number  and  static  pressures  with  COIL  gas  flow  rates.  The  design  for  a  single  injector  is 
also  discussed. 


3.2  Gas  Delivery  System 

This  experiment  uses  COIL  gas  flow  ratios  to  accurately  represent  the  fluid  mechanics  in 
a  COIL  nozzle.  In  a  COIL  nozzle,  typical  molar  gas  flows  are  4:1  helium  to  oxygen  in  the 
primary  flow,  with  a  plenum  (stagnation)  pressure  of  65  Torrf 8  Injected  with  helium  as  a  carrier 
gas,  the  iodine  molar  flow  rate  is  1.5%  of  the  oxygen  flow  rate.  The  primary  flows  are  delivered 
at  room  temperature  (300K).  To  determine  the  molar  flow  rates  of  the  primary  flow,  Eqn.  A.  10 
(from  Appendix  A)  is  used.  The  throat  and  span  of  the  nozzle  are  1  cm  and  5  cm  respectively, 
giving  a  throat  area  of  5  cm'  (see  Section  3.3  for  design  details).  The  molecular  weight  for  the 
helium  and  oxygen  mixture  is  9.6  kg/kmol  and  the  ratio  of  specific  heats  is  1.587,  giving  a  total 
molar  flow  rate  of  about  640  mmol/s,  with  514  mmol/s  of  helium  and  129  mmol/s  of  oxygen. 
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The  injected  flow  is  129  mmol/s  of  helium  and  1.5  mmol/s  of  iodine.  The  molar  flow  rate 
measurement  technique  and  the  uncertainties  are  detailed  in  Appendix  A. 

The  oxygen  and  helium  gases  flow  through  an  in-line  mixer  upstream  of  the  nozzle, 
sufficient  to  fully  mix  the  two  gases.  To  ensure  the  unifonnity  of  the  primary  gas  flow  and  to 
dissipate  any  vortical  structures  in  it,  a  flow  straightener  was  placed  downstream  of  the  in-line 
mixer  and  just  upstream  of  nozzle  plenum.  The  flow  straightener  is  a  series  of  honeycomb  and 
fine  mesh  screens  to  take  out  both  axial  and  lateral  turbulence.  The  flow  straightener  design  is 
discussed  in  Appendix  B.  The  nozzle  is  pumped  with  a  vacuum  system  that  has  a  capacity  of 
793  m  / min  and  can  pump  to  a  pressure  as  low  as  0.2  Torr. 


3.3  Nozzle  Design 

The  nozzle  was  designed  with  a  Mach  number  of  2.2  at  the  nozzle  exit  plane  (NEP). 
Mach  2.2  decreases  the  static  temperature  from  300  K  at  the  plenum  to  less  than  150  K.  The 
contour  of  the  nozzle  was  designed  using  a  computer  program  developed  at  the  University  of 

TO 

Illinois  at  Urbana-Champaign.  The  nozzle  is  contoured  on  the  top  and  bottom,  while  the 
sidewalls  are  flat  polycarbonate  plates.  The  contour  has  been  adjusted  to  compensate  for  a 
boundary  layer  build-up  along  the  contours  and  sidewalls.  A  thorough  description  of  the  nozzle 
design  and  the  boundary  layer  correction  is  given  in  Appendix  C.  The  throat  of  the  nozzle  is  1 
cm,  while  the  span  is  5  cm.  The  NEP  is  4.3  cm  downstream  of  the  nozzle  throat.  The  contour 
profile  is  then  extended  to  22.5  cm  downstream  of  the  throat  at  a  constant  2-degree  angle, 
matching  the  contour  angle  at  the  NEP.  The  nozzle  profile  from  the  plenum  to  the  exit  of  the 
nozzle  is  displayed  in  Figure  3. 
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Figure  3.  The  nozzle  contour,  from  the  plenum  to  the  exit  of  the  nozzle. 

The  nozzle  consists  of  a  stainless  steel  main  body  frame,  top  and  bottom  stainless  steel 
contour  inserts,  and  polycarbonate  sidewalls.  Multiple  inserts,  each  with  an  injector  plenum, 
were  made  so  that  different  iodine  injector  scenarios  could  be  tested.  The  injectors  were 
fabricated  by  Wire-Tech,  Inc.  after  the  inserts  were  complete,  using  sinker  EDM  (Electrical 
Discharge  Machining).  Three  holes  penetrate  the  inserts  in  the  spanwise  direction  around  the 
injector  plenum  (which  is  near  the  throat  of  the  nozzle).  Heaters  are  placed  in  these  holes  to 
keep  the  nozzle  block  warm  so  that  the  gaseous  iodine  does  not  freeze  out  in  the  plenum  or 
injectors.  Figure  4  is  a  picture  of  the  nozzle  frame  with  only  one  insert  installed,  while  Figure  5 
is  a  picture  of  the  frame  and  inserts  after  installation.  Figure  6  is  a  picture  of  the  nozzle  on  the 
test  stand,  where  the  polycarbonate  side  plates  and  pressure  taps  are  visible. 
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Figure  4.  Nozzle  frame  with  top  contour  insert  installed. 


Figure  5.  Nozzle  frame  with  top  and  bottom  contour  inserts,  shown  after  installation. 


Figure  6.  Nozzle  installed  on  test  stand.  The  polycarbonate  side  plates  are  visible  along  with  the  pressure 
taps  on  the  top  and  bottom  of  the  nozzle. 
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3.3.1  Pressure  and  Mach  Number  Measurements 


To  verify  the  nozzle  design,  pressure  measurements  were  made  without  iodine 
injection.  One  of  the  smooth  polycarbonate  walls  was  replaced  with  a  polycarbonate  wall 
that  had  fittings  at  2.5-cm  increments.  A  Pitot  tube  was  inserted  in  these  fittings.  Because 
the  flow  is  supersonic,  the  stagnation  pressures  upstream  and  downstream  of  the  shock  wave 
formed  by  the  Pitot  tube  were  used  to  calculate  the  Mach  number  in  the  flow.  The  isentropic 
relationship  for  pressure  and  Mach  number  and  the  relationship  for  Mach  numbers  up-  and 
downstream  of  a  normal  shock  were  used  to  detennine  the  Mach  number  from  the  ratio  of 


stagnation  pressures.  From  Lifshitz  and  Landau,  the  relationship  between  the  Mach  number 
upstream  of  the  normal  shock  wave,  My,  and  the  ratio  of  stagnation  pressures,  where  Poi  is 
the  stagnation  pressure  upstream  of  the  normal  shock  and  P02  is  the  stagnation  pressure 
downstream  of  the  normal  shock,  is40 


The  stagnation  pressure  upstream  of  the  nonnal  shock  was  measured  in  the  plenum  of 
the  nozzle,  whereas  the  stagnation  pressure  downstream  of  the  normal  shock  wave  was 
measured  by  the  Pitot  tube.  The  Pitot  tube  was  placed  along  the  centerline  of  the  nozzle. 

The  ratio  of  specific  heats,  y,  is  assumed  constant  (1.587)  in  this  calculation.  The  Mach 
number  distribution,  as  determined  from  the  measured  stagnation  pressures,  is  displayed  in 
Figure  7.  At  the  NEP,  the  Mach  number  is  close  to  the  design  number  of  2.2.  Along  with  the 
experimental  results,  Figure  7  displays  the  Mach  number  from  a  three  dimensional 
computational  fluid  dynamics  (CFD)  calculation  for  this  nozzle,  which  agrees  with  the 
measured  values.41 
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Distance  from  Nozzle  Throat  (mm) 
Figure  7.  Mach  number,  measured  and  calculated  with  3-D  CFD. 


Static  pressure  measurements  were  made  along  the  top  and  bottom  walls  of  the 
nozzle.  As  can  be  seen  in  Figure  6,  there  are  pressure  taps  along  the  top  and  bottom 
contours.  The  taps  are  staggered  from  top  to  bottom  and  are  2  cm  apart,  allowing  static 
pressure  measurements  to  be  made  in  1-cm  increments.  The  measured  static  pressures,  along 
with  3-D  CFD  results  and  the  isentropic  prediction  are  shown  in  Figure  8. 41  The  3-D  CFD 
calculation  matches  well  with  the  measured  results. 
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Figure  8.  Wall  static  pressures,  as  measured,  calculated  with  3-D  CFD,  and  calculated  using 
isentropic  relationships. 


3.3.2  Injector  Design 

The  contour  inserts  were  designed  so  that  injectors  were  machined  after  the  inserts 
were  manufactured.  The  first  injection  scenario  is  a  single  injector  positioned  on  the 
centerline  of  the  bottom  contour,  0.25  cm  from  the  nozzle  throat.  It  is  angled  45°  to  the 
primary  flow.  The  throat  of  the  injector  is  0.5 1  mm  in  diameter,  while  the  minor  axis  of  the 
jet  exit  (which  is  elliptical  because  of  the  wall  curvature)  is  1.0  mm.  The  area  ratio  of  the 
injector  gives  a  jet  exit  Mach  number  of  3.4.  The  45°  angle  and  the  high  Mach  number  of 
3.4  are  designed  to  provide  substantial  jet  penetration  into  the  supersonic  crossflow.  Though 
these  may  not  be  the  optimum  parameters,  this  is  the  design  to  be  studied  here.  Figure  9  is  a 
schematic  of  the  injector. 
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Figure  9.  Schematic  of  a  single  injector,  units  are  given  in  millimeters. 
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4.0  PLIF  DIAGNOSTICS 


4.1  Introduction 

The  PLIF  system  consists  of  a  Neodymium:Yttrium-Aluminum-Garnet  (Nd:YAG)  laser, 
tunable  dye  laser,  laser-sheet-forming  optics,  and  an  intensified  charge-coupled  device  (ICCD) 
camera.  This  chapter  will  present  the  laser  specifications  and  the  optical  setup  in  detail. 


4.2  Lasers 

For  this  experiment,  a  Nd:YAG  laser  is  used  to  pump  a  tunable  dye  laser.  The  Nd:YAG 
laser  is  operated  at  the  second  harmonic  (532  mn),  with  400  mJ/pulse.  It  has  a  10-Hertz  pulse 
rate  and  an  8-  to  10-ns  pulse  width.  Rhodamine  6G  is  used  in  the  dye  laser,  giving  a  laser 
wavelength  range  of  559  to  576  mn.  The  pulse  width  is  10  ns  and  the  pulse  energy  is  about  100 
mJ. 

Temperature  variation  in  the  nozzle  can  affect  the  iodine  spectral  lines  that  are  being 
pumped  and  the  iodine  fluorescence  lifetime,  both  of  which  can  impact  the  fluorescence 
intensity.  Fluorescence  intensity,  for  the  purposes  of  this  study,  is  directly  related  to  the 
concentration  of  the  injected  flow.  Because  the  temperature  variation  in  nozzle  is  considerable, 
by  as  much  as  300  K  from  the  nozzle  throat  to  the  exit,  it  is  essential  to  determine  the  wavelength 
and  bandwidth  at  which  the  laser  is  operated  and  the  fluorescence  lifetime  throughout  the  nozzle. 
To  verify  the  laser  wavelength  and  bandwidth,  the  laser  was  calibrated  with  a  Hollow  Cathode 
Discharge  Lamp  (HCL),  as  discussed  in  Section  4.2.1.  Section  4.2.2  discusses  the  method  to 
measure  the  iodine  fluorescence  lifetime  and  maps  the  fluorescence  lifetime  throughout  the 
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nozzle.  During  the  PLIF  experiments,  the  dye  laser  was  operated  at  565.095  nm  with  a 
bandwidth  of  0.024  nm.  At  this  wavelength  and  bandwidth,  eight  iodine  spectral  lines  were 
pumped.  Two  primary  spectral  lines  are  being  pumped  throughout  the  temperature  range. 
Because  the  same  lines  are  being  pumped  throughout  the  jet  of  the  nozzle,  intensity  adjustments 
are  not  required  to  accurately  compare  PLIF  images.  The  fluorescence  lifetime  has  the  largest 
variance  at  the  exit  of  the  jet,  where  the  temperature  gradient  is  the  greatest.  This  large  variance 
in  lifetime  must  be  considered  when  analyzing  the  PLIF  images. 

4.2.1  Dye  Laser  Calibration 

The  optogalvanic  effect  was  exploited  using  an  HCL  to  detennine  the  bandwidth  and 
approximate  wavelength  correction  of  the  dye  laser.  Optogalvanic  wavelength  calibration 
with  an  HCL  is  discussed  by  Zhu  et  al.  "  The  optogalvanic  effect  is  the  impedance  change  of 
a  discharge  tube  when  illuminated  by  a  second  discharge  tube  or  tunable  dye  laser.  The 
voltage  across  the  discharge  tube  is  measured  while  the  dye  laser  is  tuned  to  the  wavelength 
of  transition  of  a  species  present  in  the  discharge  of  the  lamp.  In  this  experiment,  a 
neon/gold  lamp  is  used  as  the  HCL.  Unfortunately,  neon  does  not  have  a  strong  transition 
near  the  peak  dye  laser  wavelength  of  565  nm.  However,  neon  does  have  a  strong  transition 
with  an  excitation  wavelength  of  588.1895  nm,  which  is  still  achievable  with  the  Rhodamine 
6G  dye.  The  laser  bandwidth  and  an  approximate  wavelength  correction  will  be  detennined 
at  this  wavelength.  To  detennine  the  dye  laser  wavelength  correction  at  the  laser  reading  of 
564.7077  nm  (as  used  in  the  experiments),  iodine  fluorescence  response,  as  measured  with  a 
photodiode,  is  plotted  versus  laser  reading  wavelength.  This  plot  is  compared  to  and  adjusted 
to  match  the  iodine  spectrum. 
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A  schematic  of  the  lamp  setup  is  shown  in  Figure  10.  The  lamp  is  powered  by  a  low- 
amperage,  high-voltage  power  supply.  About  350  volts  and  1  mA  were  needed  to  start  the 
lamp,  though  only  170  volts  were  needed  to  operate  the  lamp  after  startup. 

Dye  Laser 
Beam 


R  =  12.5  kQ  R  =  5.6  kQ  PS 

Figure  10.  A  schematic  of  the  hollow  cathode  lamp  (HCL)  setup.  PS  is  the  power  supply,  R  is  a 
resistor,  C  is  the  capacitor,  BC  is  the  box  controller,  and  PC  is  the  computer. 

To  capture  the  588. 1895  ntn  transition,  the  dye  laser  was  set  up  to  step  over  a 
wavelength  range  that  would  include  the  transition.  This  first  pass  was  crude  and  was  used 
to  find  a  range  of  wavelengths  that  included  the  588.1895  ntn  transition.  The  lamp  had  an 
increase  in  voltage  between  the  dye  laser  settings  of  587.9  mn  and  588.05  nm.  The  dye  laser 
was  then  setup  to  scan  from  587.9000  nm  to  588.0500  nm  in  0.0001-nm  increments.  A 
computer  recorded  the  wavelength  and  the  lamp  voltage  at  each  increment.  The  lamp 

voltage  versus  wavelength  is  plotted  in  Figure  11.  The  lamp  voltage  peaked  when  the  dye 
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laser  was  tuned  to  587.963 1  nm,  rather  than  the  published  588.1895-nm  peak  for  neon.  The 
dye  laser  wavelength,  therefore,  is  incorrect  by  0.2264  nm  at  588. 1895  nm.  The  correction 
for  the  laser  is  not  linear,  so  this  is  only  an  estimate  for  the  calibration  at  the  operating  laser 
wavelength  during  testing  (where  the  laser  is  set  at  564.7077  nm).  From  Fig.  3.2,  the 
bandwidth  is  noted  by  the  full-width,  half-maximum  (FWHM)  of  the  lamp  voltage  response. 
The  bandwidth  is  0.024  nm,  compared  to  the  specified  bandwidth  of  0.002  nm  at  570  nm.4' 
This  is  a  large  enough  bandwidth  to  excite  multiple  iodine  transitions. 


587.88  587.90  587.92  587.94  587.96  587.98  588.00  588.02  588.04  588.06 

Dye  Laser  Wavelength  (nm) 

Figure  11.  Normalized  neon  lamp  voltage  versus  dye  laser  wavelength. 


The  dye  laser  wavelength  must  be  calibrated  over  the  wavelength  used  during  the 

PLIF  tests.  Because  the  lamp  does  not  have  strong  transitions  over  the  operating  wavelength 

of  the  laser,  a  glass  cell  with  iodine  crystals  (under  vacuum)  was  set  up  to  monitor  iodine 

excitation.  The  dye  laser  passed  through  the  cell,  while  a  photodiode  (with  a  band  pass  filter 

fitted  to  only  pass  the  iodine  fluorescence)  monitored  the  iodine  fluorescence.  The  dye  laser 

wavelength  was  scanned  in  0.0002-nm  increments  over  a  large  wavelength  range.  The 
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normalized  photodiode  response  versus  dye  laser  wavelength  (after  a  0.387  nm  correction)  is 
plotted  in  Figure  12.  The  dye  laser  needed  to  be  adjusted  by  0.387  nm  to  match  the  iodine 
spectrum.  Therefore,  the  wavelength  of  564.7077  nm  used  during  the  PLIF  tests  is  actually 
565.0947  nm. 


Dye  Laser  Wavelength,  Translated  +0.387nm 

Figure  12.  Normalized  photodiode  response  versus  an  adjusted  dye  laser  wavelength.  The  normalized 
iodine  spectrum  is  also  plotted  for  comparison. 

To  identify  which  iodine  lines  are  being  excited,  the  photodiode  response  and  iodine 

spectrum  are  plotted  for  a  smaller  wavelength  range  centered  on  565.095  nm.  Figure  13 

shows  the  photodiode  response,  the  iodine  spectrum,  and  the  peaks  of  the  individual  iodine 

spectral  lines.  The  bandwidth  is  noted  so  that  the  spectral  lines  that  are  excited  are  evident. 

Within  the  bandwidth,  there  are  eight  strong  spectral  lines.  These  spectral  lines  are  defined 

in  Table  1.  The  ground  and  excited  states  for  vibration  and  rotation  are  given,  where  v"  is  the 

ground-state  (X  electronic  state)  vibration  level  and  v '  is  the  excited-state  (B  electronic  state) 
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Normalized  Value 


vibration  level.  Likewise,  J'  is  the  ground-state  rotational  level  and  /'is  the  excited  state 


rotational  level. 


Dye  Laser  Wavelength  (nm) 

Figure  13.  Photodiode  response,  iodine  spectrum,  and  individual  spectral  lines  versus  wavelength. 
The  dye  laser  bandwidth  and  operating  wavelength  are  shown. 
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Table  1.  Wavelength,  relative  absorption  values,  and  vibration  and  rotational  transitions  for  the  spectral 
lines  included  in  the  bandwidth  centered  around  565.095  nm.44 


X  (nm) 

Normalized 

X  State 

B  State 

Absorption 

v" 

J" 

v' 

J' 

565.0830 

0.07 

1 

115 

22 

116 

565.0840 

0.81 

0 

29 

18 

28 

565.0873 

0.31 

0 

129 

20 

130 

565.0918 

0.96 

1 

79 

21 

80 

565.0971 

0.35 

0 

125 

20 

124 

565.0974 

0.66 

0 

34 

18 

35 

565.0975 

1.00 

1 

75 

21 

74 

565.0977 

0.52 

0 

96 

19 

97 

The  injected  flow  temperature  ranges  from  100  K  to  400  K.  The  population  of  the 
molecules  in  a  particular  ground  state  change  with  temperature.  The  ground  state  population, 

N,  is  45 


Noz(2J  +  l)e~{G+F)/kT  (4) 

The  population  is  dependent  on  the  degeneracy  of  the  rotational  level,  (2  J+  1),  the 
energy  of  the  v  vibration  level,  G,  and  the  energy  of  the  J  rotational  level,  F.  For  a  vibrating 
rotator,  the  rovibrational  energy  is  45 
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(5) 


In  this  equation,  <x>e  is  the  equilibrium  vibration  frequency,  coexe  is  an  anharmonicity 
constant,  Bv  is  the  rotational  constant,  Dv  is  the  centrifugal  distortion  constant,  and  HOT  are 


higher-order  terms.  The  rotational  and  centrifugal  distortion  constants  are  determined  by 


45 


30 


(8) 

(9) 


The  constants  coe,  (OeXe,  Be,  and  ae  are  tabulated  in  Herzberg  for  the  iodine  molecule  in 
ground  state  as  214.57  cm'1,  0.6127  cm'1,  0.03735  cm'1,  and  0.000117  cm'1,  respectively.45 

With  the  vibration  and  rotational  levels  known,  the  population  as  a  function  of 
temperature  can  be  plotted.  Figure  14  shows  the  populations  for  the  ground-state  levels  that 
are  pumped  within  the  bandwidth  of  the  dye  laser.  Though  the  population  units  are  arbitrary, 
the  plot  illustrates  the  population  trend  with  increasing  temperature.  It  is  evident  that  the  v"  = 
0,  J'  =  29  and  v”  =  0,  J”  =  34  are  the  strongest  transition  lines  being  pumped  by  the  dye  laser 
within  the  entire  temperature  range  of  the  jet.  Figure  15  displays  the  population  fraction  as  a 
function  of  the  ground-state  rotational  level,  for  v"=0  and  varying  temperatures.  For  the 
dominant  transition  lines  pumped  by  the  laser  (J"  =  29  and  34),  the  population  fraction  can 
change  by  as  much  as  0.01  from  100K  to  500K,  which  would  not  drastically  alter  the  iodine 
fluorescence  intensity.  From  the  Boltzmann  distributions,  the  iodine  fluorescence 
wavelength  is  consistent  throughout  the  jet  of  the  nozzle  and  intensity  adjustments  are  not 
required  to  accurately  compare  PLIF  images. 
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Figure  14.  Boltzmann  distribution  for  the  ground-state  vibration  and  rotational  levels  that  are  excited 
in  the  bandwidth  of  the  dye  laser  at  565.095nm. 
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Figure  15.  Boltzmann  distribution  with  constant  temperature  and  v"=0. 
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4.2.2  Iodine  Fluorescence  Lifetime 


Another  factor  that  influences  the  collected  fluorescence  intensity  is  the  iodine 
fluorescence  lifetime.  Iodine  fluorescence  lifetime  is  most  affected  by  the  temperature  and 
pressure  in  the  nozzle,  both  of  which  influence  the  number  of  molecular  collisions.  The 
temperature  in  the  jet  of  the  nozzle  decreases  by  as  much  as  300  K  from  the  jet  exit  to  the 
nozzle  exit,  whereas  the  pressure  decreases  by  as  much  as  3500  Pa  (see  Fig.  C.5).  Capelle 
and  Broida  published  fluorescence  lifetimes  and  quenching  cross  sections  of  an  excited  state 
of  the  iodine  molecule.46  They  produced  the  following  relationship  between  the  observed 
fluorescence  lifetime  r0bs  and  pressure  and  temperature: 


1  1 

Tn 


m,  kT 


PU  +  17  He 


VnekT 


PHe  +  <*0, 


MoM 


(10) 


where  r0  is  the  natural  lifetime  of  molecular  iodine.  The  collisional  cross  sections  for  self- 


and  foreign  gas  quenching  (helium  and  oxygen)  are  represented  by  as,  ciHe,  and  ooi, 
respectively.  The  iodine,  helium,  and  oxygen  pressures  are  Pn,  PHe,  and  P02,  respectively. 
mn  is  the  mass  of  an  iodine  molecule,  p  is  the  reduced  mass  for  an  I2  foreign  gas  molecular 
collision  ( juHe  and  /U02  for  helium  and  oxygen,  respectively),  and  k  is  the  Boltzmann  constant. 
To  isolate  the  effects  that  temperature  and  pressure  have  on  the  fluorescence  lifetime,  it  is 
useful  to  use  the  ideal  gas  law,  P=pRT.  Replacing  temperature  with  PI pR  and  assuming  the 

density  is  constant  gives  a  P  relationship  to  1/ z0bs,  which  indicates  that  a  decrease  in 

1/2 

pressure  increases  the  fluorescence  lifetime.  Replacing  pressure  with  pRT  gives  a  T 
relationship  to  1/ z0bs,  which  indicates  that  a  decrease  in  temperature  would  increase  the 
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iodine  fluorescence  lifetime.  With  the  dramatic  decreases  in  pressure  and  temperature  from 
the  jet  exit  to  the  nozzle  exit,  the  fluorescence  lifetime  is  expected  to  increase. 

PLIF  imaging  was  used  to  determine  the  iodine  fluorescence  lifetime  in  different 
sections  of  the  nozzle.  The  ICCD  camera  has  the  capability  of  capturing  a  series  of  images 
with  increasing  delay  from  an  external  trigger.  Gaseous  iodine  was  seeded  in  the  primary 
flow  of  the  nozzle,  with  the  primary  gas  flow  rates  matching  normal  test  conditions.  The 
nozzle  contour  inserts  were  blank  so  that  there  were  no  injectors  present  in  this  study. 

The  camera  was  aimed  perpendicular  to  the  flow  direction,  capturing  streamwise 
images  of  the  flow  with  a  2-ns  gate.  The  camera  collected  images  at  the  same  rate  as  the 
lasers  (10  Hz)  and  was  triggered  by  the  Nd:YAG  laser.  The  camera  delay  was  set  such  that 
the  first  image  was  captured  before  the  laser  pulse.  Each  progressive  image  was  delayed  by 
1  ns,  such  that  the  final  image  had  a  total  delay  of  at  a  least  100  ns  (depending  on  the 
fluorescence  lifetime  at  that  location).  The  images  captured  the  excitation  laser,  the  most 
intense  iodine  fluorescence,  the  progressive  decrease  in  fluorescence  intensity,  and  zero 
intensity  (after  the  lifetime  of  the  iodine  fluorescence).  Figure  16  presents  the  normalized 
intensity  versus  time  for  a  single  image  pixel.  The  pixel  location  is  0.48  cm  downstream  of 
the  nozzle  throat  and  at  the  nozzle  half-height.  The  fluorescence  lifetime  is  the  time  between 
the  peak  intensity  and  1/e  fraction  of  the  peak  intensity.  At  this  location,  the  fluorescence 
lifetime  is  approximately  50  ns. 
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Figure  16.  Intensity  versus  time  for  a  single  image  pixel,  located  0.48  cm  downstream  of  the  nozzle 
throat,  at  the  nozzle  half-height.  The  blue  line  is  the  normalized  intensity  while  the  green  line  is  an 
exponential  fit  to  the  data.  The  lifetime  is  measured  from  the  peak  intensity  to  1/e  of  the  peak 
intensity. 


Figure  17  presents  the  normalized  intensity  versus  time  for  an  image  pixel  located 
5.86  cm  downstream  of  the  throat  (at  the  nozzle  half-height).  The  fluorescence  lifetime  at 
this  location  is  about  120  ns.  This  is  a  much  greater  lifetime  than  at  0.48  cm  downstream  of 
the  nozzle  throat.  To  map  fluorescence  lifetimes  versus  downstream  distance  from  the  throat 
of  the  nozzle,  the  pixel  lifetimes  were  averaged  over  the  center  of  the  nozzle  (30  pixels  along 
the  centerline  for  the  512  pixels  in  image  width)  for  each  image.  The  average  fluorescence 
lifetime  for  each  image  is  displayed  in  Figure  18.  The  fluorescence  lifetime  increases  from 
about  55  ns  to  140  ns,  which  is  significant. 
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Figure  17.  Intensity  versus  time  for  a  single  image  pixel,  located  5.86  cm  downstream  of  the  nozzle 
throat,  at  the  nozzle  half-height.  The  blue  line  is  the  normalized  intensity  while  the  green  line  is  an 
exponential  fit  to  the  data.  The  lifetime  is  measured  from  the  peak  intensity  to  1/e  of  the  peak 
intensity. 
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Figure  18.  Two-point  moving  average  of  the  iodine  fluorescence  lifetime  versus  downstream  distance 
from  throat.  The  lifetime  is  the  average  for  an  image  field  of  view  along  the  centerline,  where  the 
distance  from  the  nozzle  throat  is  measured  at  the  image  center. 


To  assess  the  collected  intensity  of  an  image,  the  area  under  the  exponential  curve  for 
a  20-ns  gate,  10  ns  after  the  peak  intensity,  can  be  calculated.  Using  intensities  that  have  not 
been  normalized  by  the  maximum  image  intensity,  the  areas  are  33,285  ns-count  for  the  0.48- 
cm  location  and  67,257  ns-count  at  the  5.86-cm  location.  From  the  normalized  images, 
where  2775  counts  is  the  maximum  intensity  for  the  0.48-cm  location  and  4668  counts  for 
the  5.86-cm  location,  the  areas  are  12  ns-normalized  counts  and  14.4  ns-normalized  counts, 
respectively.  Figure  19  displays  the  integrated  normalized  lifetime  versus  distance  from  the 
nozzle  throat.  Each  data  point  is  averaged  along  the  centerline  of  each  image,  as  was  done  in 
Figure  18.  However,  an  error  bar  is  shown  to  depict  plus/minus  one  standard  deviation  from 
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the  image  mean.  From  this  plot,  there  is  not  a  strong  correlation  between  the  integrated 
normalized  lifetime  and  downstream  distance  from  the  nozzle  throat. 
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Figure  19.  Averaged  integrated  normalized  fluorescence  lifetime  versus  downstream  distance  from 
the  throat.  The  error  bars  are  plus/minus  one  standard  deviation  of  the  integrated  normalized  lifetime 
for  each  image. 


In  collecting  the  lifetime  images,  there  were  fluctuations  in  the  molar  flow  rate  of  the 
iodine  in  the  flow.  This  fluctuation  could  affect  the  lifetime  data,  so  the  integrated 
nonnalized  lifetime  was  plotted  versus  the  average  iodine  molar  flow  rate,  shown  in  Figure 
20.  There  is  not  a  strong  correlation  between  the  integrated  normalized  lifetime  and  iodine 
molar  flow  rate. 
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Figure  20.  Averaged  integrated  normalized  fluorescence  lifetime  versus  iodine  molar  flow  rate. 


The  lifetime  data  presented  thus  far  is  with  the  iodine  seeded  in  the  primary  flow  of 
the  nozzle.  The  iodine  that  is  present  in  the  jet  flow  is  injected  at  a  higher  temperature  of  430 
K  and  static  pressure  of  about  1600  Pa.  The  jet  flow  also  has  a  greater  concentration  of 
iodine,  which  could  result  in  greater  self-quenching.  To  accurately  represent  the 
fluorescence  lifetime  in  the  jet,  lifetime  measurements  were  made  with  a  single  injector  and 
iodine  seeded  only  in  the  injector.  A  contour  plot  of  fluorescence  lifetime  for  four  different 
downstream  locations  is  shown  in  Figure  21 . 
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Figure  21.  Contour  plots  of  the  fluorescence  lifetimes  (color  bar  units  are  nanoseconds)  at  four 
overlapping  locations  in  the  nozzle.  Spatial  dimension  units  in  plots  are  in  mm,  x-location  is  measured 
from  the  nozzle  throat,  and  v-location  is  measured  from  the  nozzle  half-height. 


Fluorescence  lifetime  was  obtained  over  the  entire  length  of  the  nozzle.  Close  to  the 
injector,  the  location  and  height  of  the  jet  do  not  fluctuate  during  data  collection,  but  these 
begin  to  fluctuate  with  increasing  downstream  distance  due  to  inherent  unsteadiness  of  the  jet 
flow.  Because  of  this  variance  in  location  and  height,  the  fluorescence  lifetime  data  is 
distorted.  Figure  21  gives  a  good  representation  of  the  fluorescence  lifetime  range  in  the 
injected  flow.  The  largest  lifetime  gradient  is  seen  at  the  exit  of  the  injector,  where  the 
lifetime  ranges  from  50  to  170  ns.  This  variation  is  most  likely  from  the  large  temperature 
variation  at  the  jet  exit.  A  3-D  computational  fluid  dynamic  (CFD)  simulation  of  this  nozzle 
has  been  done  by  Madden  et  ah,  where  the  temperature  variation  is  plotted  within  this  single 
jet  (shown  in  Figure  22).47  The  temperature  is  the  lowest  in  the  core  of  the  jet  but  increases 
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toward  the  outer  edges.  As  the  jet  issues  downstream,  the  temperature  remains  relatively 
constant,  where  the  fluorescence  lifetime  remains  within  the  range  of  150  to  170  ns. 


Figure  22.  3-D  Navier-Stokes  simulation  results  for  temperature  contours  within  a  2-D  planar  slice 
through  the  center  of  the  jet.47 


As  with  the  fluorescence  lifetime  calculations  for  the  primary  flow  seeded  with 
iodine,  plotting  the  integrated  normalized  lifetime  will  give  a  better  picture  of  what  the  ICCD 
camera  actually  sees.  Figure  23  presents  the  integrated  normalized  lifetimes  at  the  same  four 
downstream  locations  as  Figure  21.  The  integrated  normalized  lifetime  is  relatively  constant 
in  the  core  of  the  jet,  but  is  significantly  less  at  the  outer  edges  of  the  jet  near  the  orifice. 

This  variance  in  the  integrated  normalized  lifetime  is  most  likely  caused  by  the  increased 
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temperatures  (and  decreased  fluorescence  lifetimes)  at  those  edges.  This  variation  at  the  jet 
exit  should  be  considered  when  analyzing  the  PLIF  images. 
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Figure  23.  Contour  plots  of  the  integrated  normalized  lifetimes  (color  bar  units  are 
nanoseconds*normalized  lifetime)  at  four  overlapping  locations  in  the  nozzle.  Spatial  dimension  units 
in  plots  are  in  mm,  x- location  is  measured  from  the  nozzle  throat,  and  v-location  is  measured  from  the 
nozzle  half-height. 


4.3  Optics 

The  beam  of  the  NdiYAG  laser  is  aimed  directly  into  the  dye  laser.  The  beam  from  the 
dye  laser  is  turned,  using  a  BK7  prism,  into  a  series  of  three  lenses  to  collimate  and  expand  the 
beam  into  a  rectangular  sheet  with  minimal  thickness  (~600  pm).  A  schematic  of  the  laser  and 
optics  setup  is  shown  in  Figure  24,  where  the  lasers  and  camera  are  positioned  to  acquire  oblique 
images  in  the  flow.  Spherical  lenses  are  used  to  collimate  the  beam;  each  has  a  focal  length  of  8 
cm.  The  collimated  beam  then  passes  through  a  cylindrical  lens,  with  a  focal  length  of  14  cm,  to 


42 


form  a  laser  sheet  with  minimal  thickness  and  a  height  of  ~4  cm.  Two  mirrors  are  used  to  direct 
the  beam  into  the  nozzle  at  a  45-degree  angle.  Before  the  laser  sheet  enters  the  nozzle,  it  passes 
through  another  cylindrical  lens,  with  a  focal  length  of  30  cm,  to  further  decrease  the  sheet 
thickness. 
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Figure  24.  A  schematic  of  the  laser  and  optics  setup.  The  lenses  are  annotated  with  SL  for  spherical  lens, 
CL  for  cylindrical  lens,  and  M  for  mirror.  In  this  schematic,  the  lasers  and  optics  are  situated  to  acquire 
images  in  the  oblique  view. 

Along  with  oblique  images,  streamwise  images  were  captured.  Figure  25  is  a  schematic 
of  the  laser  and  optics  setup  for  streamwise  imaging.  The  laser  enters  through  a  window  in  the 
vacuum  line  and  is  aimed  down  the  centerline  of  the  nozzle,  in  the  opposite  direction  of  the  gas 


flow.  The  laser  sheet  is  -600  pm  thick  and  is  the  height  of  the  nozzle.  The  camera  is  positioned 
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perpendicular  to  the  laser  sheet  and  is  on  a  translation  stage  to  easily  obtain  pictures  from  the 
injector  to  the  nozzle  exit. 
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Figure  25.  A  schematic  of  the  lasers  and  optics  setup  to  acquire  streamwise  images.  The  laser  sheet  is 
aimed  up  the  back  end  of  the  nozzle  through  a  window. 


The  Princeton  Instruments  Acton  ICCD  camera  has  a  Generation  III  photo  intensifies 
with  GaAs  as  the  photocathode  material.4s  The  camera  is  able  to  gate  as  fast  as  2  ns  and  can 
operate  at  a  10-Hz  frequency.  The  camera  resolution  is  5 12  x  5 12  pixels.  During  the  tests,  a 
Sigma  50  mm  F2.8  Macro  camera  lens  is  used,  allowing  a  resolution  of  -0.03  mm/pixel  for 
streamwise  images  and  -0.05  mm/pixel  for  oblique  images.  The  camera  gate  time  was  set  to  20 
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ns,  which  is  sufficient  to  collect  enough  fluorescence  yet  is  fast  enough  to  capture  minute  gas 
movement.  At  the  exit  of  the  nozzle,  where  the  gases  are  traveling  with  the  highest  velocity  of 


approximately  900m/s  (see  Fig.  C.4),  an  element  in  the  gas  translates  only  0.018  mm  in  20  ns. 
This  movement  is  less  than  a  pixel  width. 
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5.0  PLIF  IMAGE  ANALYSIS  OF  SINGLE  INJECTOR 


5.1  Introduction 

The  first  nozzle  design  studied  is  a  nozzle  with  a  single  injector  on  the  bottom  wall  of  the 
nozzle  (as  described  in  Chapter  3).  To  visualize  the  flow,  PLIF  images  were  collected  in 
ensembles  of  100.  The  images  were  taken  at  progressive  downstream  locations  with  two 
different  camera  views  (streamwise  or  oblique).  The  laser  and  camera  configurations  for  both 
the  oblique  and  streamwise  views  are  shown  in  Figures  24  and  25,  respectively.  The  streamwise 
view  is  a  side  view  of  the  jet  flowing  from  left  to  right,  where  images  were  taken  along  the 
nozzle  centerline.  The  oblique-view  images  the  jet  at  a  45-degree  angle,  displaying  the  spanwise 
structure  of  the  jet. 

The  injected  flow  is  highlighted  by  the  iodine  fluorescence.  The  flow  with  the  highest 
fluorescence  intensity  is  the  unmixed  counter-rotating  vortex  pair,  while  the  less  intense  portion 
of  the  flow  is  a  mixed  portion  of  the  primary  and  injected  gasses,  and  the  portion  of  the  image 
with  no  intensity  is  the  primary  flow. 

Ensemble  averages  were  computed  for  both  the  streamwise  and  oblique  views. 

Histogram  analysis  was  performed  on  the  ensemble  averages  of  both  views  to  map  the  intensity 
variance,  which  lends  insight  into  the  mixing  quality  between  the  injected  and  primary  flows. 
With  the  oblique  images,  the  structure  size  is  measured  to  quantify  the  jet  growth  and  diffusive 
properties. 


46 


5.2  Ensemble  Averaging 


Averages  were  taken  of  each  image  ensemble  (of  100  images)  for  both  the  streamwise 
and  oblique  PLIF  images. 

5.2.1  Streamwise  Images 

A  schematic  of  the  streamwise  locations  of  the  PLIF  images  is  displayed  in  Figure  26, 
and  ensemble  averages  of  the  streamwise  images  are  displayed  in  Figures  27  through  31.  In 
each  image,  the  flow  is  from  left  to  right,  the  nozzle  wall  is  drawn  in  red,  the  x-value  is  the 
distance  from  the  nozzle  throat,  and  the  v-value  is  the  distance  from  the  nozzle  half-height. 
The  images  display  the  single  injector  from  the  nozzle  throat  to  9.6  cm  downstream  of  the 
throat.  Intensity  streaks  in  the  images  are  caused  by  optical  aberrations  when  forming  the 
excitation  laser-sheet.  Most  of  the  streamwise  images  were  taken  with  a  camera  gain  of  200, 
though  the  streamwise  image  in  Figure  29  was  taken  with  a  camera  gain  of  100.  To  account 
for  this  difference,  the  overlapping  section  of  the  images  in  Figures  28  and  29  (30  -  36  mm) 
was  used  to  determine  an  intensity  multiplier.  Each  ensemble  average  was  normalized  by  the 
maximum  intensity  at  the  exit  of  the  injector  (24,072  counts). 

In  general,  with  progressive  downstream  distance,  the  jet  is  penetrating  farther  into 
the  centerline  of  the  nozzle.  After  the  nozzle  exit  plane  (NEP),  the  jet  began  to  fluctuate 
vertically  in  time.  This  fluctuation  can  be  seen  in  a  sequence  of  instantaneous  images  at  a 
streamwise  location,  where  the  jet  height  is  continuously  changing.  The  fluctuations  become 
more  intense  with  downstream  distance.  Ensemble  averaging  the  fluctuating  images  does  not 
accurately  convey  the  maximum  jet  penetration.  Figure  30  is  evidence  of  the  jet  fluctuation 
where  the  penetration  is  much  less  than  in  Figures  29  and  3 1 . 
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Figure  27.  Ensemble  average  of  the  streamwise  images  at  the  jet  exit.  Note  that  the  iodine  fluorescence  is 
reflected  onto  the  bottom  wall  of  the  nozzle,  shown  under  the  red  line.  Color  bar  units  are  pixel  intensity. 
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Figure  28.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~22  mm  to  ~36  mm  downstream  of 
the  nozzle  throat.  Color  bar  units  are  pixel  intensity. 
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Figure  29.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~29  mm  to  ~43  mm  downstream  of 
the  nozzle  throat.  Color  bar  units  are  pixel  intensity. 
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Figure  30.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~53  mm  to  ~67  mm  downstream  of 
the  nozzle  throat.  Color  bar  units  are  pixel  intensity. 
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Figure  31.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~82  mm  to  ~96  mm  downstream  of 
the  nozzle  throat.  Color  bar  units  are  pixel  intensity. 


5.2.2  Oblique  Images 

To  study  the  spanwise  structure  of  the  jet,  oblique-view  images  were  taken.  Figure 
32  is  a  schematic  of  the  axes  used  in  notating  the  oblique  view.  The  schematic  is  a  top-view 
of  the  nozzle  with  gas  flowing  from  top  to  bottom.  The  throat,  shown  with  a  dashed  line,  is 
where  x  =  0.  Thc  v-axis  is  out-of  the  page,  starting  at  the  nozzle  half-height.  The  x45- axis 
begins  at  the  x-axis  and  follows  the  laser  sheet.  The  camera  view  is  typically  24  to  30  mm 
wide,  centered  at  the  intersection  of  the  x45-  and  x-axes.  Because  of  the  oblique  view,  the 
left-hand-side  of  each  image  is  downstream  of  the  right-hand- side.  If  the  structure  size  is  4 
mm  wide,  the  left-hand-side  is  2.8  nun  downstream  of  the  right-hand-side. 
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Figure  32.  Schematic  of  the  x-,  x45-,  and  _v-axcs  with  respect  to  the  nozzle  and  camera.  The  ,v45-axis  is 
along  the  direction  of  the  laser  sheet. 


Ensemble  averages  of  the  oblique  images  are  displayed  in  Figures  33  through  39. 

The  upstream  to  downstream  separation  of  each  45°  plane  is  approximately  10  mm,  from  13 
mm  to  75  mm  downstream  from  the  nozzle  throat.  The  oblique  images  were  taken  with  a 
camera  gain  of  100.  The  ensemble  averages  are  normalized  by  the  maximum  intensity  in 
Figure  33  (5898  counts).  The  top  and  bottom  walls  of  the  nozzle  are  indicated  by  red  dashes. 
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Figure  33.  Ensemble  average  of  the  oblique  images  at  v~13  mm.  Color  bar  units  are  pixel  intensity. 
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Figure  34.  Ensemble  average  of  the  oblique  images  at  x~23  mm.  Color  bar  units  are  pixel  intensity. 
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Figure  35.  Ensemble  average  of  the  oblique  images  at  x~35  mm.  Color  bar  units  are  pixel  intensity. 
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Figure  36.  Ensemble  average  of  the  oblique  images  at  x~45  mm.  Color  bar  units  are  pixel  intensity. 
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Figure  37.  Ensemble  average  of  the  oblique  images  at  x~55  mm.  Color  bar  units  are  pixel  intensity. 
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Figure  38.  Ensemble  average  of  the  oblique  images  at  x~65  mm.  Color  bar  units  are  pixel  intensity. 
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Figure  39.  Ensemble  average  of  the  oblique  images  at  x~15  mm.  Color  bar  units  are  pixel  intensity. 

A  pair  of  counter-rotating  vortices  is  evident  in  the  oblique  images.  In  Figures  34 
through  Figure  39,  a  recirculation  region  below  the  counter-rotating  vortex  pair  (CRVP)  is 
visible  as  well.  The  height  of  the  CRVP  is  increasing  with  downstream  distance,  along  with 
the  recirculation  region  of  the  jet.  From  Figure  39,  it  is  evident  that  the  jet  structure  extends 
to  the  centerline  of  the  nozzle  by  x  =  75  mm. 


5.3  Histograms 

To  quantify  the  mixing  between  the  jet  and  the  primary  flow,  a  histogram  analysis  was 
performed  on  both  the  streamwise  and  oblique  images.  Histograms  display  the  probability 
density  of  the  fluorescence  intensity.  For  an  unmixed  flow,  a  histogram  will  typically  be 
bimodal.  It  will  have  a  zero-intensity  peak  indicating  the  primary,  non-fluorescing  flow  and  a 
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high-intensity  peak  indicating  the  injected,  fluorescing  flow.  There  also  may  be  a  third 
intermediate-intensity  peak  if  the  primary  and  secondary  flows  are  mixing. 

The  intensities  of  the  peaks  will  shift  with  downstream  distance.  As  mixing  occurs,  the 
intermediate-  and  high-intensity  peaks  decrease  with  downstream  distance.  This  indicates  that 
the  primary  flow  is  being  advected  into  the  CRVP  where  it  is  mixed  with  the  jet,  thus  reducing 
the  fluorescence  intensity  of  the  injected  flow. 

5.3.1  Streamwise  Images 

Histograms  for  the  ensemble-averaged  streamwise  images  are  shown  in  Figure  40. 

The  normalized  intensity  is  displayed  on  the  abscissa,  while  the  number  of  pixels  at  a  specific 
intensity  is  displayed  on  the  ordinate.  The  downstream  location  of  each  histogram  is 
displayed  on  the  x-axis,  where  the  downstream  location  of  the  streamwise  image  is  taken  as 
the  location  of  the  center  pixel.  In  all  the  histograms,  a  zero-intensity  peak  is  present, 
representing  the  primary  flow.  The  high-intensity  peak  in  each  histogram,  representing 
iodine,  occurs  at  lower  intensity  levels  with  downstream  distance.  A  third,  intermediate- 
intensity  peak  fonns  at  x  =  29  mm  (red  line)  and  moves  to  lower  intensity  levels  with 
downstream  distance.  Because  of  the  decrease  in  intensity  levels  with  the  high- intensity  and 
intennediate-intensity  peaks,  mixing  is  evident  with  the  streamwise  PLIF  images. 
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Normalized  Intensity 

Figure  40.  Histograms  of  the  ensemble-averaged  streamwise  images  from  Figures  27  to  3 1 .  The  intensity 
of  each  image  is  normalized  by  the  maximum  intensity  in  Figure  27.  The  x-location  is  measured  at  the 
center  of  the  image. 


5.3.2  Oblique  Images 

As  with  the  streamwise  images,  a  histogram  analysis  was  performed  on  the  ensemble- 
averaged  oblique  images.  Histograms  for  Figures  33  through  39  are  displayed  in  Figure  41. 
The  zero-intensity  and  high-intensity  peaks  are  discernible,  but  an  intermediate-intensity 
peak  is  difficult  to  distinguish.  The  high-intensity  peaks  are  decreasing  with  downstream 
distance,  as  with  the  streamwise  images,  indicating  mixing  between  the  primary  and  injected 
flows. 
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Normalized  Intensity 

Figure  41.  Histograms  of  the  ensemble-averaged  oblique  images  from  Figures  33  to  39.  The  intensity  of 
each  image  is  normalized  by  the  maximum  intensity  in  Figure  33. 


5.4  Structure  Size 

If  diffusive  mixing  between  the  primary  and  injected  flows  occurs  with  downstream 
distance,  the  jet  structure  will  grow  in  size.  This  size  can  be  measured  by  its  area  or  length  of 
interface  between  the  jet  structure  and  primary  flow.49 

To  determine  the  area  of  the  oblique  structure,  the  grayscale  PLIF  image  is  converted  to 
monochrome,  with  a  threshold  of  30%  of  the  maximum  image  intensity.  Figure  42  is  a 
monochrome  rendition  of  Figure  39  (x  =  75  mm).  The  black  pixels,  representing  all  intensities 
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greater  than  or  equal  to  30%  of  the  maximum  image  intensity,  are  counted,  where  the  count  is 
converted  to  square  millimeters  to  give  the  structure  area.  The  area  was  determined  for  x 
positions:  13  mm,  23  mm,  35  mm,  45  mm,  55  mm,  65  mm,  and  75  mm.  Structure  area  versus 
downstream  distance  from  the  nozzle  throat  is  displayed  in  Figure  43.  The  area  increases  with 
downstream  distance,  where  the  area  more  than  doubles  from  the  x  =  13  mm  position  to  the  x  = 
75  mm  position. 
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Figure  42.  Black  and  white  image  of  the  oblique  structure  at  x~75  mm. 
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Figure  43.  Areas  of  the  oblique  image  structures  versus  downstream  distance  from  the  nozzle  throat. 

Outlines  from  the  monochrome  images  were  used  to  estimate  the  interface  length  of  the 
structures  in  the  oblique  view.  These  were  computed  using  the  edge-detection  function  Outline 
in  the  public-domain  imaging  software,  IrnageJ.50  This  function  examines  each  pixel  in  the 
monochrome  image,  and  flips  it  from  black  to  white  if  all  pixels  surrounding  it  are  black.  This 
leaves  black  pixels  only  at  the  edge  of  the  structure,  defining  its  outline.  Figure  44  is  an  outline 
of  the  structure  shown  in  Figure  39.  The  number  of  black  pixels  in  the  outline  was  counted  and 
then  converted  to  interface  length  (in  millimeters)  using  the  pixel  resolution  (0.05  mm/pixel). 
Interface  length  versus  downstream  distance  from  the  nozzle  throat  is  displayed  in  Figure  45.  As 
with  structure  area,  interface  length  increases  with  downstream  distance. 
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Figure  44.  Outline  of  the  black  and  white  image  from  Figure  42. 
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Figure  45.  interface  lengths  of  the  oblique  image  structures  versus  downstream  distance  from  the  nozzle 
throat. 
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With  both  interface  length  and  structure  area  increasing  with  downstream  distance,  it  is 
evident  that  the  CRVP  is  increasing  in  size  with  downstream  distance.  This  indicates  that 
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diffusive  mixing  between  the  injected  and  primary  flows  of  the  nozzle  is  occurring.  It  is  also 
suggestive  that  there  is  advection  of  the  main  flow  into  the  CRVP. 

The  next  chapter  will  discuss  the  investigation  of  other  injection  scenarios,  where  trip  jets 
are  used  to  perturb  the  flow.  Along  with  histograms,  interface  length  and  structure  areas  will  be 
used  to  compare  the  relative  mixing  for  each  of  the  injection  scenarios. 
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6.0  MULTIPLE  INJECTORS 


6.1  Introduction 

To  enhance  mixing  between  the  primary  and  injected  flows,  small  injectors  were 
positioned  in  the  wake  of  the  single  injector.  These  small  injectors  are  termed  trip  jets  because 
they  are  intended  to  disturb,  or  trip,  the  counter-rotating  vortex  pair  of  the  single  injector.  The 
trip  jets  are  0.25  mm  in  diameter,  angled  perpendicular  to  the  nozzle  wall,  and  have  a  jet  exit 
Mach  number  of  1 .  The  diameter  of  each  trip  jet  was  designed  to  be  significantly  smaller  than 
the  injector,  which  has  an  effective  diameter  of  1.02  mm.  With  the  sizeable  thickness  of  the 
nozzle  between  the  wall  and  injector  plenum,  the  smallest  diameter  that  could  be  manufactured 
by  Wire-Tech,  Inc.,  was  0.25  mm.  Therefore,  the  trip  jet  diameters  are  0.25  mm.  The  constant 
diameter  of  the  trip  jets,  from  the  plenum  to  the  jet  exit,  creates  the  Mach  1  design.  They  are 
angled  perpendicular  to  the  nozzle  wall,  or  83°  to  the  primary  flow,  to  allow  penetration  to  the 
injector  jet. 

The  trip  jets  were  placed  in  the  wake  of  the  injector,  where  the  flow  is  subsonic,  to 
facilitate  the  trip  jet  penetration  to  the  counter-rotating  vortex  pair  (CRVP).  The  first  injection 
scenario  tested  (Trip  Scenario  A)  consists  of  a  single  injector  placed  along  the  centerline  of  the 
nozzle,  with  its  leading  edge  1.00  mm  downstream  of  the  trailing  edge  of  the  injector,  which  is 
one  injector  diameter.  This  places  the  center  of  the  trip  jet  4.41  mm  downstream  of  the  nozzle 
throat.  Once  this  scenario  was  tested,  Trip  Scenario  B  was  created  to  determine  if  decreasing  the 
distance  between  the  injector  and  trip  would  have  an  increased  impact  on  the  disturbance  of  the 
CRVP.  The  distance  between  the  trailing  edge  of  the  injector  and  leading  edge  of  the  trip  was 
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approximately  halved  to  0.45  mm,  where  the  center  of  the  trip  jet  is  3.96  mm  from  the  nozzle 
throat.  Once  Trip  Scenario  B  was  tested,  Trip  Scenario  C  was  created  to  determine  if  two 
injectors,  3.96  mm  from  the  nozzle  throat,  each  0.33  mm  from  the  centerline,  could  destabilize 
the  CRVP  or,  at  least,  push  the  counter-rotating  vortices  apart.  The  centerline  of  the  trip  jets  was 
positioned  0.33  mm  from  the  centerline  to  hit  the  jet  on  the  outer  edge  of  each  vortex  core.  From 
Figure  33,  the  outer  edge  of  each  vortex  core  (though  it  appears  as  one  vortex  close  to  the  jet 
exit)  is  approximately  one-third  of  the  jet  width  (0.33  mm)  from  the  centerline.  Figures  46 
through  48  are  schematics  of  the  injector  and  trip  jets  from  Trip  Scenarios  A  through  C.  Table  2 
summarizes  the  dimensions  of  the  Trip  Scenarios. 


Figure  46.  Top-  and  side-view  schematic  of  the  injector  and  trip  jet  for  Trip  Scenario  A.  Units  are  in 
millimeters  and  the  primary  flow  is  left  to  right. 
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Figure  47.  Top-  and  side-view  schematic  of  the  injector  and  trip  jet  for  Trip  Scenario  B.  Units  are  in 
millimeters  and  the  primary  flow  is  left  to  right. 


Figure  48.  Top-  and  side-view  schematic  of  the  injector  and  trip  jet  for  Trip  Scenario  C.  Units  are  in 
millimeters  and  the  primary  flow  is  left  to  right. 
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Table  2.  Trip  Scenario  specifications. 


Trip  Scenario  A 

Trip  Scenario  B 

Trip  Scenario  C 

Number  of  Trip  Jets 

1 

1 

2 

Trip  Jet  Diameter 

0.25  mm 

0.25  mm 

0.25  mm 

Trip  Jet  Mach  Number 

1 

1 

1 

Distance  from 
Nozzle  Throat 

4.41  mm 

3.96  mm 

3.96  mm 

Edge-to-Edge  Distance 
from  Injector 

1.00  mm 

0.45  mm 

0.45  mm 

Distance  from 
Nozzle  Centerline 

0.00  mm 

0.00  mm 

0.33  mm 

To  compare  the  nozzles  with  the  trip  jets  to  the  nozzle  with  the  single  injector, 
streamwise  and  oblique  PLIF  images  were  taken.  Histogram  analyses  were  perfonned  on  the 
ensemble-averaged  streamwise  and  oblique  images  for  each  trip  jet  scenario.  Jet  penetration  was 
measured  with  the  ensemble-averaged  streamwise  images  while  the  interface  length  and 
structure  areas  of  the  oblique  structures  were  determined  for  comparison  to  the  nozzle  with  the 
single  injector. 

6.2  Streamwise-View  Ensemble  Averages 

As  in  Chap.  5,  streamwise  images  were  collected  in  ensembles  of  100.  For  each  trip  jet 
scenario,  the  streamwise  PLIF  images  are  normalized  by  the  maximum  intensity  at  the  exit  of  the 
injector  for  that  scenario.  The  streamwise  images  overlap  by  approximately  7  mm  and  are  taken 
from  the  jet  exit  to  approximately  75  mm  downstream  of  the  nozzle  throat.  The  PLIF  images  for 
Trip  Scenarios  A  and  B  were  taken  with  a  camera  gain  of  200,  while  the  images  for  Trip 
Scenario  C  were  taken  with  a  camera  gain  of  100.  Because  the  images  are  normalized,  they  can 
be  directly  compared.  In  each  image,  the  flow  is  from  left  to  right,  the  nozzle  wall  is  drawn  in 
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red,  the  x-value  is  the  distance  from  the  nozzle  throat,  and  the  v-value  is  the  distance  from  the 
nozzle  half-height. 


6.2.1  Trip  Scenario  A 

Ensemble-averaged  streamwise  PLIF  images  for  Trip  Scenario  A  (shown  in  Figure 
46)  are  displayed  in  Figures  49  through  58.  The  camera  gain  is  200  for  all  images  and  each 
image  is  normalized  by  the  maximum  intensity  at  the  injector  exit.  Figure  49  is  a  mosaic  of 
all  streamwise  images,  excluding  Figure  54  (where  the  jet  height  did  not  match  well  with  the 
upstream  and  downstream  images). 
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Figure  49.  Ensemble  average  of  the  streamwise  images  at  the  jet  exit  for  Trip  Scenario  A.  Note  that  the 
iodine  fluorescence  is  reflected  from  the  bottom  wall  of  the  nozzle,  shown  under  the  red  line.  Color  bar 
units  are  pixel  intensity. 
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Figure  50.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~8  mm  to  ~22  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  A.  Color  bar  units  are  pixel  intensity. 

14000 

12000 

10000 

8000 

6000 

4000 

2000 


16  17  18  19  20  21  22  23  24  25  26  27  28 
x  (mm) 

Figure  51.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~15  mm  to  ~29  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  A.  Color  bar  units  are  pixel  intensity. 
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Figure  52.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~22  mm  to  ~36  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  A.  Color  bar  units  are  pixel  intensity. 
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Figure  53.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~29  mm  to  ~43  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  A.  Color  bar  units  are  pixel  intensity. 
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Figure  54.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~36  mm  to  ~50  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  A.  Color  bar  units  are  pixel  intensity. 
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Figure  55.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~43  mm  to  ~57  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  A.  Color  bar  units  are  pixel  intensity. 
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Figure  56.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~49  mm  to  ~63  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  A.  Nozzle  wall  is  not  present  in  image.  Color  bar  units  are  pixel 
intensity. 
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Figure  57.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~57  mm  to  ~7 1  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  A.  Nozzle  wall  is  not  present  in  image.  Color  bar  units  are  pixel 
intensity. 
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Figure  58.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~64  mm  to  ~78  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  A.  Nozzle  wall  is  not  present  in  image.  Color  bar  units  are  pixel 
intensity. 
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6.2.2  Trip  Scenario  B 


Ensemble-averaged  streamwise  PLIF  images  for  Trip  Scenario  B  (illustrated  in 
Figure  47)  are  displayed  in  Figures  60  through  69.  The  camera  gain  is  200  for  all  images  and 
each  image  is  normalized  by  the  maximum  intensity  at  the  injector  exit.  Figure  70  is  a 
mosaic  of  all  streamwise  images. 


-3  -2  -1  0  1  2  3  4  5  6  7  8  9  10  11 


4 

x  10 


5 

4 

3 

2 

1 

0 


x  (mm) 

Figure  60.  Ensemble  average  of  the  streamwise  images  at  the  jet  exit,  for  Trip  Scenario  B.  Color  bar  units 
are  pixel  intensity. 
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Figure  61.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~4  mm  to  ~18  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  B.  Color  bar  units  are  pixel  intensity. 
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Figure  62.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~1 1  mm  to  ~25  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  B.  Color  bar  units  are  pixel  intensity. 
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Figure  63.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~18  mm  to  ~32  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  B.  Color  bar  units  are  pixel  intensity. 
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Figure  64.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~25  mm  to  ~39  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  B.  Color  bar  units  are  pixel  intensity. 
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Figure  65.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~32  mm  to  ~46  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  B.  Color  bar  units  are  pixel  intensity. 
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Figure  66.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~39  mm  to  ~53  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  B.  Color  bar  units  are  pixel  intensity. 
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Figure  67.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~45  mm  to  ~60  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  B.  Color  bar  units  are  pixel  intensity. 
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Figure  68.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~53  mm  to  ~67  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  B.  Color  bar  units  are  pixel  intensity. 


79 


60  61  62  63  64  65  66  67  68  69  70  71  72  73  74 


5000 

4000 

3000 

2000 

1000 

0 


(mm) 

Figure  69.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~60  mm  to  ~74  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  B.  Color  bar  units  are  pixel  intensity. 
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6.2.3  Trip  Scenario  C 


Ensemble-averaged  streamwise  PLIF  images  for  Trip  Scenario  C  (illustrated  in 
Figure  48)  are  displayed  in  Figures  71  through  81.  The  camera  gain  is  100  for  all  images  and 
each  image  is  normalized  by  the  maximum  intensity  at  the  injector  exit.  Figure  82  is  a 
mosaic  of  all  streamwise  images. 

16000 
14000 
12000 
10000 
8000 
6000 
4000 
2000 


-3  -2  -1  0  1  2  3  4  5  6  7 
x  (mm) 

Figure  71.  Ensemble  average  of  the  streamwise  images 
are  pixel  intensity. 
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Figure  72.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~5  mm  to  ~19  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  C.  Color  bar  units  are  pixel  intensity. 
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Figure  73.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~9  mm  to  ~24  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  C.  Color  bar  units  are  pixel  intensity. 
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Figure  74.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~16  mm  to  ~3 1  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  C.  Color  bar  units  are  pixel  intensity. 
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Figure  75.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~22  mm  to  ~37  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  C.  Color  bar  units  are  pixel  intensity. 
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Figure  76.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~30  mm  to  ~45  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  C.  Color  bar  units  are  pixel  intensity. 
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Figure  77.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~36  mm  to  ~52  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  C.  Color  bar  units  are  pixel  intensity. 
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Figure  78.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~42  mm  to  ~58  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  C.  Color  bar  units  are  pixel  intensity. 
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Figure  79.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~49  mm  to  ~65  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  C.  Color  bar  units  are  pixel  intensity. 
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Figure  80.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~56  mm  to  ~72  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  C.  Color  bar  units  are  pixel  intensity. 
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Figure  81.  Ensemble  average  of  the  streamwise  images  of  the  jet  from  ~62  mm  to  ~77  mm  downstream  of 
the  nozzle  throat,  for  Trip  Scenario  C.  Color  bar  units  are  pixel  intensity. 
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6.3  Oblique- View  Ensemble  Averages 


For  each  trip  jet  scenario,  oblique-view  PLIF  images  were  taken  as  described  in  Section 
5.2.2.  The  images  were  taken  in  ensembles  of  100  and  averaged.  For  each  trip  jet  scenario,  the 
oblique  images  were  normalized  by  the  maximum  intensity  in  the  most-upstream  oblique  image 
(closest  to  the  jet  exit).  When  present,  the  top  and  bottom  walls  of  the  nozzle  are  indicated  by 
red  dashes. 


6.3.1  Trip  Scenario  A 

Ensemble-averaged  oblique  PLIF  images  for  Trip  Scenario  A  (see  Figure  46)  are 
displayed  in  Figures  83  through  90.  The  camera  gain  is  100  for  all  images  except  Figures  85 
and  86,  where  the  intensities  are  adjusted  as  was  done  in  Figure  29.  Each  image  is 
normalized  by  the  maximum  intensity  in  Figure  83. 
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Figure  83.  Ensemble  average  of  oblique  images  at  ,r~10  mm,  for  Trip  Scenario  A.  Color  bar  units  are 
pixel  intensity. 
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Figure  84.  Ensemble  average  of  oblique  images  at  x~i 3  mm,  for  Trip  Scenario  A.  Color  bar  units  are 
pixel  intensity. 
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Figure  85.  Ensemble  average  of  oblique  images  at  x~25  mm,  for  Trip  Scenario  A.  Color  bar  units  are 
pixel  intensity. 
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Figure  86.  Ensemble  average  of  oblique  images  atx~35  mm,  for  Trip  Scenario  A.  Color  bar  units  are 
pixel  intensity. 
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Figure  87.  Ensemble  average  of  oblique  images  at  x~45  mm,  for  Trip  Scenario  A.  Color  bar  units  are 
pixel  intensity. 
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Figure  88.  Ensemble  average  of  oblique  images  at  x~55  mm,  for  Trip  Scenario  A.  Color  bar  units  are 
pixel  intensity. 
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Figure  89.  Ensemble  average  of  oblique  images  at  x~65  mm,  for  Trip  Scenario  A.  Color  bar  units  are 
pixel  intensity. 
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Figure  90.  Ensemble  average  of  oblique  images  at  x~l 5  mm,  for  Trip  Scenario  A.  Color  bar  units  are 
pixel  intensity. 


6.3.2  Trip  Scenario  B 


Ensemble-averaged  oblique  PLIF  images  for  Trip  Scenario  B  (see  Figure  47)  are 
displayed  in  Figure  91  through  97.  The  camera  gain  is  100  for  all  images  except  Figures  96 
and  97,  where  the  intensities  are  adjusted  as  was  done  in  Figure  29.  Each  image  is 
normalized  by  the  maximum  intensity  in  Figure  91. 

This  trip  jet  design  produced  some  unsteadiness  in  the  jet.  To  ensemble-average  the 
oblique  structure,  the  structure  in  each  image  was  centered  using  a  method  described  by 
Vorobieff  et  al.51  For  each  instantaneous  image,  the  center  of  intensity  (x/,  v/)  is  determined 
in  the  same  fashion  the  center  of  gravity  would  be  calculated  in  solid  mechanics: 
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where  xtj  is  the  x-pixel,  ytj  is  the  v-pixel,  and  Iq  is  the  intensity  at  pixel  (i,j).  Each 
instantaneous  image  was  then  translated  to  the  center  pixel  for  each  of  the  100  images  in  the 
ensemble.  The  ensemble  was  then  averaged.  The  centering  allows  for  a  better  representation 
of  the  oblique  structure  with  fluctuating  instantaneous  images.  The  fluctuating  ensembles 
were  centered,  but  the  steady  ensembles  were  not. 
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Figure  91.  Ensemble  average  of  oblique  images  at  jc~1  2.3  mm,  for  Trip  Scenario  B.  Color  bar  units  are 
pixel  intensity. 
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Figure  92.  Ensemble  average  of  oblique  images  at  x~25  mm,  for  Trip  Scenario  B.  Instantaneous  images 
were  centered.  Color  bar  units  are  pixel  intensity. 
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Figure  93.  Ensemble  average  of  oblique  images  of  the  jet  at  *~35  mm,  for  Trip  Scenario  B.  Instantaneous 
images  were  centered.  Color  bar  units  are  pixel  intensity. 
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Figure  94.  Ensemble  average  of  oblique  images  at  x~45  mm,  for  Trip  Scenario  B.  Instantaneous  images 
were  centered.  Color  bar  units  are  pixel  intensity. 
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Figure  95.  Ensemble  average  of  oblique  images  at  x~55  mm,  for  Trip  Scenario  B.  Instantaneous  images 
were  centered.  Color  bar  units  are  pixel  intensity. 
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Figure  96.  Ensemble  average  of  oblique  images  of  the  jet  at  x~65  mm,  for  Trip  Scenario  B.  Instantaneous 
images  were  centered.  Color  bar  units  are  pixel  intensity. 
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Figure  97.  Ensemble  average  of  oblique  images  at  x~l 5  mm,  for  Trip  Scenario  B.  Instantaneous  images 
were  centered.  Color  bar  units  are  pixel  intensity. 
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6.3.3  Trip  Scenario  C 


Ensemble-averaged  oblique  PLIF  images  for  Trip  Scenario  C  (see  Figure  48)  are 
displayed  in  Figures  98  through  105.  The  camera  gain  is  100  for  all  images  and  each  image 
is  normalized  by  the  maximum  intensity  in  Figure  98.  Trip  Scenario  C  also  produced  some 
unsteadiness  in  the  jet.  The  ensembles  with  unsteady  instantaneous  images  were  centered,  as 
described  in  Section  6.3.2. 

10000 
9000 
8000 
7000 
6000 
5000 
4000 
3000 
2000 
1000 
0 

-12-10  -8  -6  -4  -2  0  2  4  6  8  10 

x45  (mm) 

Figure  98.  Ensemble  average  of  oblique  images  at  x~l  mm,  for  Trip  Scenario  C.  Color  bar  units  are  pixel 
intensity. 
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Figure  99.  Ensemble  average  of  oblique  images  at  x~10.5  mm,  for  Trip  Scenario  C.  Color  bar  units  are 
pixel  intensity. 
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Figure  100.  Ensemble  average  of  oblique  images  at  x~25  mm,  for  Trip  Scenario  C.  Instantaneous  images 
were  centered.  Color  bar  units  are  pixel  intensity. 
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Figure  101.  Ensemble  average  of  oblique  images  at  x~35  mm,  for  Trip  Scenario  C.  Instantaneous  images 
were  centered.  Color  bar  units  are  pixel  intensity. 
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Figure  102.  Ensemble  average  of  oblique  images  at  x~45  mm,  for  Trip  Scenario  C.  instantaneous  images 
were  centered.  Color  bar  units  are  pixel  intensity. 
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Figure  103.  Ensemble  average  of  oblique  images  at  x~55  mm,  for  Trip  Scenario  C.  Color  bar  units  are 
pixel  intensity. 
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Figure  104.  Ensemble  average  of  oblique  images  at  x~65  mm,  for  Trip  Scenario  C.  Color  bar  units  are 
pixel  intensity. 
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Figure  105.  Ensemble  average  of  oblique  images  at  x~75  mm,  for  Trip  Scenario  C.  Instantaneous  images 
were  centered.  Color  bar  units  are  pixel  intensity. 


6.4  Streamwise-View  Histograms 

A  histogram  analysis  was  performed  on  the  ensemble-averaged  streamwise  images  for 
Trip  Scenarios  A,  B,  and  C.  The  analysis  is  described  in  detail  in  Section  5.3.  The  normalized 
intensity  is  displayed  on  the  x-axis,  while  the  number  of  pixels  at  a  specific  intensity  is  displayed 
on  the  z-axis.  The  downstream  location  of  each  histogram  is  displayed  on  the  y-axis,  where  the 
downstream  location  of  the  streamwise  image  is  taken  as  the  location  of  the  center  pixel. 

6.4.1  Trip  Scenario  A 

Histograms  for  the  ensemble-averaged  streamwise  images  for  Trip  Scenario  A  are 
shown  in  Figure  106.  In  all  the  histograms,  a  zero-intensity  peak  is  present,  representing  the 
primary  flow.  The  high-intensity  peak  (representative  of  the  injected  flow)  in  each  histogram 
is  decreasing  in  intensity  with  downstream  distance.  The  intensity  decrease  is  greater  than 
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the  decrease  shown  in  the  histograms  for  the  single  injector  (shown  in  Figure  40).  This 
greater  intensity  decrease  signifies  that  the  trip  jet  enables  greater  mixing  between  the 
primary  and  injected  flows. 


Figure  106.  Histograms  of  the  ensemble-averaged  streamwise  images  from  Figures  49  to  58,  Trip 
Scenario  A.  The  intensity  of  each  image  is  normalized  by  the  maximum  intensity  in  Figure  49.  The  x- 
location  is  measured  at  the  center  of  the  image. 


6.4.2  Trip  Scenario  B 

Histograms  for  the  ensemble-averaged  streamwise  images  for  Trip  Scenario  B  are 
displayed  in  Figure  107.  The  second,  high-intensity  peak  has  an  even  greater  decrease  in 
intensity  than  the  high-intensity  peaks  from  Figure  106,  which  is  evidence  that  Trip  Scenario 
B  has  even  greater  mixing  than  Trip  Scenario  A. 
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Normalized  Intensity 

Figure  107.  Histograms  of  the  ensemble-averaged  streamwise  images  from  Figures  60  to  69,  Trip 
Scenario  B.  The  intensity  of  each  image  is  normalized  by  the  maximum  intensity  in  Figure  60.  The  x- 
location  is  measured  at  the  center  of  the  image. 


6.4.3  Trip  Scenario  C 

Histograms  for  the  ensemble-averaged  streamwise  images  for  Trip  Scenario  C  are 
shown  in  Figure  108.  The  high-intensity  peak  has  a  greater  decrease  in  intensity  than  with 
Trip  Scenario  A,  but  perhaps  not  as  great  of  a  decrease  as  Trip  Scenario  B.  This  signifies 
that  Trip  Scenario  C  has  greater  mixing  between  the  primary  and  injected  flows  than  Trip 
Scenario  A,  but  no  greater  than  Trip  Scenario  B. 
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Normalized  Intensity 

Figure  108.  Histograms  of  the  ensemble-averaged  streamwise  images  from  Figures  71  to  81,  Trip 
Scenario  C.  The  intensity  of  each  image  is  normalized  by  the  maximum  intensity  in  Figure  71.  The  x- 
location  is  measured  at  the  center  of  the  image. 


6.5  Oblique-View  Histograms 

A  histogram  analysis  was  performed  on  the  ensemble-averaged  oblique  images  for  the 
three  different  trip  jet  scenarios.  While  the  high-intensity  peaks  generally  decrease  in  intensity 
with  downstream  distance  for  each  trip  jet  scenario,  the  pattern  of  this  decrease  does  not  match 
that  of  the  streamwise-view  histograms.  The  oblique-view  histograms  indicate  that  Trip 
Scenario  A  has  the  greatest  mixing  between  the  primary  and  injected  flows,  whereas  the 
streamwise-view  histograms  indicate  that  Trip  Scenario  B  has  the  greatest  mixing. 
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6.5.1  Trip  Scenario  A 


Histograms  for  the  ensemble-averaged  oblique  images  for  Trip  Scenario  A  are  shown 
in  Figure  109.  The  high-intensity  peak  has  a  greater  decrease  in  intensity,  with  downstream 
distance,  than  the  nozzle  with  the  single  injector.  This  indicates  that  this  trip  jet  scenario  has 
greater  mixing  between  the  primary  and  injected  flows  than  the  nozzle  with  no  trip  jets. 


Normalized  Intensity 

Figure  109.  Histograms  of  the  ensemble-averaged  oblique  images  from  Figures  83  to  90,  Trip  Scenario  A. 
The  intensity  of  each  image  is  normalized  by  the  maximum  intensity  in  Figure  83. 

6.5.2  Trip  Scenario  B 

Histograms  for  the  ensemble-averaged  oblique  images  for  Trip  Scenario  B  are 
displayed  in  Figure  1 10.  Fromx  =  12.3  mm  to  x  =  55  mm,  the  high-intensity  peak  decreases 
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in  intensity  with  downstream  distance.  At  x  =  65  mm  and  x  =  75  mm,  the  high-intensity  peak 
actually  increases  in  intensity.  Perhaps  this  is  because  of  the  trip  jet  unsteadiness  and 
misrepresentation  of  the  jet  structure  (even  though  the  images  were  centered  before  ensemble 
averaging). 


Normalized  Intensity 

Figure  110.  Histograms  of  the  ensemble-averaged  oblique  images  from  Figures  91  to  97,  Trip  Scenario  B. 
The  intensity  of  each  image  is  normalized  by  the  maximum  intensity  in  Figure  91. 


6.5.3  Trip  Scenario  C 

Histograms  for  the  ensemble-averaged  oblique  images  for  Trip  Scenario  C  are 
displayed  in  Figure  111.  The  high-intensity  peak  decreases  in  intensity  from  x  =  7  mm  to  x  = 
65  mm,  but  increases  at  x  =  75  mm.  This  is  most  likely  from  the  jet  unsteadiness,  as  with  the 


oblique-view  histograms  for  Trip  Scenario  B. 
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Normalized  Intensity 

Figure  111.  Histograms  of  the  ensemble-averaged  oblique  images  from  Figures  98  to  105,  Trip  Scenario 
C.  The  intensity  of  each  image  is  normalized  by  the  maximum  intensity  in  Figure  98. 


6.6  Structure  Size 

The  structure  area  and  interface  lengths  of  the  oblique-view  structures,  for  Trip  Scenarios 
A,  B,  and  C  were  measured.  The  method  for  measuring  the  area  and  interface  length  of  a 
structure  is  detailed  in  Section  5.4.  The  threshold  for  converting  the  grayscale  PLIF  images  to 
monochrome  was  set  at  30%,  as  was  done  in  Chap.  5.  Structure  area  versus  downstream 
distance  is  displayed  in  Figure  112.  The  values  for  each  trip  jet  scenario  and  the  single  jet  are 
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shown  for  comparison.  Close  to  the  jet  exit  (x  <  25  mm),  the  single  jet  has  a  comparable  area  to 
the  jet  with  trip  jet(s).  Trip  Scenarios  A  and  C  have  greater  structure  areas  between  x  =  35  and  x 
=  65,  while  Trip  Scenario  B  does  not  have  as  significant  of  an  impact  on  the  structure  area.  Atx 
=  75  mm,  the  structure  areas  for  Trip  Scenarios  A  and  C  decrease  because  of  the  jet  fluctuations 
caused  by  the  trip  jets,  which  decreases  the  structure  size  in  the  ensemble-averaged  image. 
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Figure  112.  Structure  area  versus  downstream  distance  for  the  single  jet  and  Trip  Scenarios  A,  B,  and  C. 

Figure  113  displays  the  interface  length  versus  downstream  distance  for  the  single 
injector  and  the  three  trip  jet  scenarios.  Close  to  the  exit  of  the  injector,  the  oblique  structure 
interface  lengths  are  similar  for  the  single  jet  and  trip  scenarios.  With  increasing  downstream 
distance,  Trip  Scenarios  A  and  C  produce  an  oblique  structure  with  much  greater  interface  length 
than  that  of  the  single  jet,  except  at  x  =  75  mm.  The  decrease  in  structure  area  and  interface 
length  at  x  =  75  mm  could  be  caused  by  the  jet  unsteadiness.  The  ensemble-averaged  oblique 
structures  for  Trip  Scenarios  B  and  C  atx  =  75  nun  (Figures  97  and  105,  respectively)  appear 
significantly  smaller  and  distorted  from  the  averages  at  upstream  positions,  even  though  the 
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instantaneous  images  have  been  centered.  This  unsteadiness  could  cause  the  ensemble  average 
(with  centering)  to  misrepresent  the  structure  area  and  interface  length  at  downstream  locations. 
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Figure  113.  Interface  length  versus  downstream  distance  for  the  single  jet  and  Trip  Scenarios  A,  B,  and  C. 
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6.7  Jet  Penetration 

With  the  ensemble-averaged  streamwise  PLIF  images,  jet  penetration  was  measured. 

The  penetration  was  measured  as  the  y-location  of  the  top  of  the  jet.  The  streamwise  images 
were  converted  to  monochrome,  with  a  threshold  of  30%  maximum  intensity.  The  top  of  the  jet, 
at  each  x-location,  was  then  determined.  For  the  single  jet  and  the  three  trip  jet  scenarios,  the  jet 
penetration  versus  downstream  distance  is  displayed  in  Figure  114.  Although  jet  penetration 
generally  increases  in  the  x-direction,  there  is  a  decrease  in  penetration  with  the  single  jet  at  60 
mm,  with  Trip  Scenario  A  at  50  mm,  with  Trip  Scenario  B  at  60  mm,  and  Trip  Scenario  C  at  65 
mm.  The  decrease  in  penetration  is  most  likely  caused  by  jet  unsteadiness  in  the  nozzle.  The 
single  jet  and  the  three  jet  scenarios  produced  a  jet  that  eventually  crossed  the  nozzle  centerline. 
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Trip  Scenario  B  crossed  the  centerline  before  the  other  cases  and  had  the  farthest  penetration  in 
the  nozzle. 
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Figure  114.  Jet  penetration  versus  downstream  distance  for  the  single  jet  and  Trip  Scenarios  A,  B,  and  C. 
Penetration  measurements  are  taken  from  the  ensemble-averaged  streamwise  images. 
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7.0  CONCLUSION  AND  RECOMMENDATIONS 


7.1  Summary  of  Results  and  Conclusions 

In  this  study,  a  supersonic  nozzle  with  COIL  applications  was  designed  and  tested  using 
PLIF  imaging.  The  nozzle  was  designed  with  a  Mach  number  of  2.2  at  the  nozzle  exit  plane 
(NEP),  where  the  downstream  distance  from  the  nozzle  throat,  x,  was  43  mm.  Pitot  tube 
measurements  and  3-D  computational  fluid  dynamic  (CFD)  computations  confirmed  a  Mach 
number  close  to  the  design  number  at  this  location. 

The  focus  of  this  study  was  to  examine  the  influence  trip  jets  may  exert  on  the  flow 
structure  and  mixing.  To  elucidate  the  trip  jet  effect  on  the  flow,  four  different  iodine  injection 
scenarios  were  imaged  with  the  PLIF  technique.  The  first  scenario  is  a  single  injector  placed 
downstream  of  the  nozzle  throat.  In  an  attempt  to  increase  mixing  between  the  injected  and 
primary  flows,  trip  jets  were  placed  downstream  of  the  single  injector  with  the  intention  to 
destabilize  the  counter-rotating  vortex  pair  (CRVP)  of  the  single  jet  or  at  least  push  the  jet  farther 
into  the  centerline  of  the  nozzle.  Three  different  trip  jet  scenarios  were  tested:  1)  Trip  Scenario 
A,  where  a  single  trip  jet  (0.25  mm  in  diameter)  is  placed  one  injector  diameter  (1.00  mm) 
downstream  of  the  single  injector,  2)  Trip  Scenario  B,  where  the  single  trip  jet  from  Trip 
Scenario  A  is  placed  closer  to  the  single  injector  (0.45  mm),  and  3)  Trip  Scenario  C,  where  two 
trip  jets  are  placed  in  the  same  downstream  location  as  Trip  Scenario  B,  but  each  are  0.33  mm 
from  the  nozzle  centerline.  Each  injection  scenario  was  imaged  in  the  streamwise  and  oblique 
views. 


112 


From  the  PLIF  images,  four  different  analyses  were  performed  on  the  images.  Histogram 
analyses  were  performed  on  both  the  streamwise-  and  oblique-view  ensemble-averaged  images. 
From  each  of  the  oblique-view  images,  the  area  of  the  oblique  structure  and  the  structure 
interface  length  were  measured.  The  final  analysis  was  the  measurement  of  the  jet  penetration 
from  the  ensemble-averaged  streamwise-view  images.  With  these  four  measurement/analysis 
techniques,  the  mixing  quality  and  penetration  can  be  compared  for  the  four  different  injection 
scenarios. 

The  histogram  analysis  was  perfonned  on  both  the  streamwise-  and  oblique-view  images. 
Histograms  plot  the  probability  density  of  the  normalized  fluorescence  intensity.  With  a 
histogram,  there  are  at  least  two  peaks,  the  zero-intensity  peak  represents  the  primary,  non¬ 
fluorescing  flow,  while  the  high-intensity  peak  represents  the  fluorescing,  injected  flow. 
Histograms  for  each  downstream  location  were  plotted  on  the  same  axis  for  each  injection 
scenario.  When  mixing  between  the  primary  and  injected  flows  occur,  the  normalized  intensity 
of  the  high-intensity  peak  decreases  with  downstream  distance  (indicating  a  dilution  of  the 
fluorescing  species).  A  greater  decrease  in  nonnalized  intensity  signifies  greater  mixing  between 
the  two  flows.  From  the  streamwise-view  images,  all  three  trip  scenarios  had  a  greater  intensity 
decrease  in  the  high-intensity  peak  than  with  the  single  injector.  With  the  single  injector,  the 
greatest  decrease  in  nonnalized  intensity  for  the  high-intensity  peak  was  from  1.0  to  0.2,  whereas 
Trip  Scenarios  A,  B,  and  C  had  peaks  that  decreased  to  0. 12,  0.06,  and  0.04,  respectively.  This 
signifies  that  Trip  Scenarios  B  and  C  had  greater  mixing  than  Trip  Scenario  A  and  the  single 
injector.  From  the  oblique-view  images,  the  high-intensity  peaks  decreased  to  0.2  for  the  single 
injector  and  0.1,  0.13,  and  0.04  for  Trip  Scenarios  A,  B,  and  C,  respectively.  From  the  oblique- 
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view  histograms,  Trip  Scenario  C  has  the  greatest  mixing,  while  Trip  Scenarios  A  and  B  had 
comparable  mixing  and  the  single  injector  had  the  least  amount  of  mixing. 

From  the  oblique-view  PLIF  images,  the  structure  size  and  length  of  interface  between 
the  injected  and  primary  flows  were  measured.  The  images  were  converted  to  monochrome, 
with  a  30%  intensity  threshold,  to  facilitate  area  and  interface  length  measurement.  For  all 
injection  scenarios,  the  structure  area  increased  with  downstream  distance.  Trip  Scenario  C  had 
the  greatest  area  of  36  mm  at  x  =  65  mm.  The  maximum  areas  for  the  single  injector,  Trip 
Scenario  A,  Trip  Scenario  B,  and  Trip  Scenario  C  were  23  mm”,  33  mm”,  19  mm  ,  and  23  mm”, 
respectively.  The  structure  areas  decreased  at  x  =  65  mm  for  Trip  Scenario  B  and  atx  =  75  mm 
for  Trip  Scenarios  A  and  C.  This  decrease  in  area  may  be  caused  by  the  jet  fluctuations.  The 
interface  length  measurement  has  the  same  trends  as  the  structure  area.  Trip  Scenarios  A  and  C 
had  the  longest  length  at  x  =  65  mm,  with  62  mm  and  69  mm,  respectively.  As  with  structure 
area,  the  interface  length  decreased  at  jc  =  65  mm  for  Trip  Scenario  B  and  at  75  mm  for  Trip 
Scenarios  A  and  C.  Because  of  the  jet  fluctuations  at  the  downstream  locations,  it  is  difficult  to 
interpret  this  data. 

With  a  larger  structure  area  and  interface  length,  the  contact  between  the  injected  and 
primary  flows  is  increased.  This  greater  contact  is  indicative  of  diffusive  mixing  in  the  nozzle. 
However,  in  some  cases,  the  jet  is  fluctuating  so  much  that  the  ensemble-averaged  structure  area 
is  smaller  than  the  structures  in  each  instantaneous  image.  The  fluctuations  dilute  the  edges  of 
the  jet  structure  in  the  ensemble  averages.  The  fluctuations  actually  enable  greater  contact 
between  the  primary  and  injected  flows,  which  increases  the  mixing  capability.  With  the  area 
and  interface  length  plots  (Figures  1 12  and  1 13),  the  measurements  are  only  useful  in  the  non¬ 
fluctuating  portions  of  the  nozzle,  up  to  55  mm.  From  the  structure  area  measurements,  the  Trip 
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Scenarios  A,  B,  and  C  all  have  better  mixing  than  the  single  injector,  with  Trip  Scenario  C 
having  the  best  results.  The  interface  length  results  also  indicate  that  the  trip  jets  increased  the 
mixing,  but  Trip  Scenario  B  has  the  best  results  (greatest  length).  These  results  are  comparable 
with  the  histogram  results,  where  Trip  Scenarios  B  and  C  had  better  mixing  than  the  single 
injector  and  Trip  Scenario  A. 

The  final  analysis  with  the  PLIF  images  is  the  jet  penetration  measurement  from  the 
streamwise  ensemble-averaged  images.  Jet  penetration  is  the  v-location  of  the  top  of  the  jet, 
which  is  defined  from  the  monochrome  image  (where  the  intensity  threshold  is  set  at  30%). 

From  the  jet  penetration  plot  (Figure  1 14),  Trip  Scenario  B  has  a  consistently  greater  penetration 
than  Trip  Scenarios  A  and  C  and  the  single  injector.  Interestingly,  the  jet  penetration  with  Trip 
Scenario  A  is  less  than  with  the  single  injector.  Jet  fluctuations  are  evident  with  Trip  Scenarios 
B  and  C  at  x  ~  57-60  mm. 

From  the  PLIF  images,  the  mixing  between  the  primary  and  injected  flows  was  measured 
by  four  different  methods:  1)  histogram  analysis,  2)  oblique  structure  area,  3)  oblique  structure 
interface  length,  and  4)  jet  penetration.  From  these  four  analyses,  it  is  evident  that  the  trip  jets 
did  influence  the  single  injector.  The  effects  of  the  trip  jet(s)  included  altered  mixing  (as 
measured  with  the  histograms),  considerable  change  in  oblique  structure  size,  and  differences  in 
the  interface  length  between  the  injected  and  primary  flows,  as  wells  as  a  noticeable  change  in 
the  jet  penetration.  The  trip  jets  also  increased  the  jet  fluctuations  at  a  downstream  distance,  x, 
greater  than  55  mm  (though  most  notable  with  Trip  Scenarios  B  and  C).  Trip  Scenarios  B  and  C 
have  an  increased  jet  penetration  and  greatest  mixing  when  measured  with  the  streamwise-view 
histograms.  Because  of  the  jet  fluctuations,  Trip  Scenario  B  has  a  smaller  oblique  structure  area, 
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oblique  structure  interface  length,  and  less  reduction  in  intensity  with  the  oblique  histograms 
than  Trip  Scenarios  A  and  C. 

In  comparing  the  results  of  the  three  trip  jet  scenarios,  it  is  evident  that  the  close 
placement  of  the  trip  jets  to  the  single  injector  in  Trip  Scenarios  B  and  C  increases  the  jet 
penetration  more  than  the  location  of  the  trip  jet  in  Trip  Scenario  A.  Mixing,  as  measured  with 
the  streamwise-view  histograms,  is  also  greatest  with  Trip  Scenarios  B  and  C.  Because  the  jet 
fluctuations  are  most  notable  with  the  oblique-view  images,  it  is  difficult  to  compare  the  area, 
interface  length,  and  oblique -view  histograms  in  detennining  the  best  trip  jet  scenario  for  the 
greatest  mixing.  With  the  increased  jet  fluctuations  with  Trip  Scenario  B,  when  compared  to 
Trip  Scenario  C,  a  single  trip  jet  just  downstream  of  the  injector  may  be  more  desirable  than  two 
offset  injectors  in  destabilizing  the  injected  flow.  These  studies  of  the  trip  jets  are  by  no  means 
exhaustive,  but  they  serve  a  dual  purpose,  first  by  experimentally  confirming  their  effectiveness, 
with  relatively  small  energy  input  into  the  flow,  and  second,  by  providing  experimental 
benchmarks  for  future  numerical  studies  of  the  problem.  Some  of  the  likely  directions  for  future 
work  are  summarized  in  the  following  section. 

7.2  Recommendations  for  Future  Work 

The  effects  of  single  or  multiple  trip  jets  on  a  single  injector  in  supersonic  flow  have  been 
investigated.  Three  trip  jet  scenarios  have  been  imaged  and  analyzed.  To  fully  understand  what 
drives  jet  unsteadiness,  jet  penetration,  and  mixing  between  the  primary  and  injected  flows,  more 
trip  jet  scenarios  need  to  be  tested.  In  this  study,  a  single  trip  jet,  placed  on  the  nozzle  centerline, 
was  tested  at  two  different  downstream  locations.  The  trip  jet  placed  closest  to  the  injector 
increased  the  jet’s  penetration  and  mixing.  Testing  more  distances  would  determine  an  optimal 
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trip  jet  distance.  Also,  the  use  of  double  trip  jets  in  Trip  Scenario  C  needs  optimization.  From 
the  oblique  PLIF  images,  it  appears  that  the  trip  jets  hit  the  outside  of  the  counter-rotating 
vortices  (see  Figure  99),  rather  than  the  center  of  each  vortex,  as  intended.  Decreasing  the 
distance  from  the  nozzle  centerline  may  produce  the  desired  unsteadiness  in  the  CRVP.  Upon 
optimization,  testing  entire  rows  of  injectors  and  trip  jets  would  reveal  the  mixing  capability  in 
the  nozzle  and  flow  unifonnity  (which  is  desirable  for  applications  to  a  COIL  system). 
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Appendix  A:  Flow  Rate  Measurement  and  Uncertainty 


A.l  Introduction 

This  appendix  details  the  molar  flow  rate  measurement  and  the  uncertainty  analysis  for 
the  measurement  of  the  primary  and  carrier  flow  rates,  including  the  flow  rate  of  iodine. 
Thermocouples,  pressure  transducers,  orifices,  and  photometers  were  all  used  in  the  flow 
measurements  and  each  have  an  uncertainty  which  will  be  included  in  this  analysis. 

The  uncertainty  of  the  iodine  molar  flow  rate  is  important  when  analyzing  the  PLIF 
images.  The  fluorescence  intensity  is  directly  affected  by  the  iodine  concentration  in  the  flow. 
Uncertainty  in  the  iodine  molar  flow  rate  or  iodine  concentration  influences  the  interpretation  of 
the  PLIF  images. 


A.2  Method  to  Measure  Molar  Flow  Rate 

To  measure  the  carrier  helium,  primary  helium,  and  primary  oxygen  molar  flow  rates, 
orifices  were  used  to  create  a  choke  point  in  their  flow.  Pressure  and  temperature  were  measured 
upstream  of  the  choke  point,  where  the  flow  is  considered  to  be  stagnate  (Mach  number  less  than 
0.05).  With  the  stagnation  pressure  and  temperature  measurements,  the  molecular  weight  of  the 
gas,  and  the  ratio  of  specific  heat  capacities  for  the  gas,  the  molar  flow  rate  can  be  determined 
from  the  isentropic  relationships,  conservation  of  mass,  and  ideal  gas  law.  An  orifice  schematic 
is  displayed  in  Figure  A.l,  below. 
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Figure  A.l.  A  schematic  of  an  orifice,  with  the  flow  from  left  to  right.  The  pressure  and  temperature 
upstream  of  the  orifice  are  stagnation  values. 


To  calculate  the  molar  flow  rate,  start  with  the  conservation  of  mass  equation,  Eqn.  A.l, 
where  p  is  the  gas  density,  A  is  the  orifice  area,  V  is  the  flow  velocity,  and  m  is  the  mass  flow 
rate. 

m  =  pAV  (A.l) 

From  the  Ideal  Gas  Law, 

P  =  pRT  (A.  2) 

where  P  is  pressure,  T  is  temperature,  and  R  is  the  gas  constant.  Equation  A.  1  can  then  be 
written  as 


PAY 


(A. 3) 


The  velocity  is  related  to  the  Mach  number,  M,  by  the  speed  of  sound,  a,  as 

V  =  Ma  (A.4) 

The  speed  of  sound  is  a  function  of  temperature  and  the  ratio  of  specific  heats  for  the  gas 

a  -  j yRT  (A. 5) 
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Equations  A.4  and  A. 5  reduce  Eqn.  A. 3  to 


til  =  PAM 


(A.  6) 


At  the  orifice,  the  Mach  number  is  one.  The  temperature  and  pressure  cannot  be 
measured  at  the  orifice,  but  could  be  deduced  from  the  stagnation  values  by  using  the  isentropic 
equations.  Also,  the  area  of  the  orifice  needs  to  be  adjusted  for  boundary  layer  build  up. 

Orifices  are  calibrated  for  the  actual  flow  rate  based  on  the  Reynolds  number.  They  are  usually 
calibrated  by  the  manufacturer  and  given  a  value  for  the  discharge  coefficient,  Cci,  where  unity  is 
an  orifice  without  a  boundary  layer. 

The  isentropic  equations  for  static  temperature  and  pressure  are 

(A.  7) 


T  =  T0  1  +  — — -M2 

L  2 


P 


(A.  8) 


With  the  isentropic  relationships,  a  Mach  number  of  one  at  the  orifice,  and  the  Cci  value 
incorporated  to  allow  for  the  reduction  in  orifice  area,  Eqn.  A. 6  becomes 


(A. 9) 


Equation  A. 9  gives  the  mass  flow  rate,  but  the  molar  flow  rate  is  desired.  The  molecular  weight, 
M,  of  the  gas  is  divided  through  to  give  molar  flow  rate.  Also,  the  gas  constant  is  replaced  by 
the  Universal  Gas  Constant,  St,  divided  by  molecular  weight.  Therefore,  the  molar  flow  rate  is 
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(A.  10) 


With  Eqn.  A.  10,  the  molar  flow  rate  of  a  flow  can  be  determined  by  the  gas  molecular  weight, 
ratio  of  specific  heats,  orifice  area  and  discharge  coefficient,  and  the  pressure  and  temperature 
upstream  of  the  orifice. 


A.3  Primary  and  Carrier  Molar  Flow  Rate  Uncertainties 

When  detennining  the  molar  flow  rate,  uncertainty  may  arise  from  the  accuracy  of  the 
measured  pressure  and  temperature  of  the  gas  and  from  the  measured  discharge  coefficient  of  the 
orifice  (which  takes  into  account  the  uncertainty  in  the  orifice  diameter).  For  simplicity  in  the 
uncertainty  calculations,  Eqn.  A.  10  can  be  reduced  to 

:V  -  fA  /y  (A.11) 

where 


(A.  12) 


The  uncertainty  in  the  molar  flow  rate  is  given  by  AN  ,  where  A  refers  to  the  uncertainty  of  a 
specified  measurement. 


(A.  13) 


where 


^L  =  I^b 

dCd  ^ 


(A.  14) 
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(A.  15) 


™=C^B 

dPo  ^ 

8N  _  1  Cd  B 
dT0  2  r03/2 


(A.  16) 


The  uncertainty  in  the  pressure  transducers  is  0.25%.  The  uncertainty  in  the 
thennocouples  is  a  conservative  1%  on  the  Kelvin  scale.  The  discharge  coefficient,  Cd,  has 
uncertainty  from  a  few  different  factors.  An  orifice  is  usually  calibrated  with  air.  Using  the 
orifice  with  another  gas  species  introduces  uncertainty  in  the  measured  discharge  coefficient 
value.  The  uncertainty  of  the  discharge  coefficient  with  the  calibrated  species  is  0.25%.  From 
Johnson  et  ah,  using  a  gas  that  differs  from  the  calibration  gas  can  change  the  discharge 
coefficient  by  2%.  From  the  calibration  uncertainty,  the  fact  that  the  gas  species  is  different 
than  that  used  for  calibration,  and  the  uncertainty  in  the  orifice  diameter,  the  orifice  discharge 
coefficient  uncertainty  was  set  at  a  conservative  5%. 

Table  A.  1  displays  a  sample  of  measured  and  calculated  values  used  in  calculating  the 
flow  rate  uncertainties  for  carrier  and  primary  helium  and  primary  oxygen.  In  each  case,  the 
flow  rate  uncertainty  is  5%,  with  the  uncertainty  in  the  orifice  discharge  coefficient  being  the 
dominant  factor. 
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Table  A.l.  Sample  calculation  of  carrier  and  primary  gas  flow  rates  and  uncertainties. 


Carrier  Helium 

Primary  Helium 

Primary  Oxygen 

Orifice  Diameter 

0.055  in 

1 .40  mm 

0.070  in 

1.78  mm 

0.075  in 

1.91  mm 

y 

1.66 

1.66 

1.40 

M 

4  g/mol 

4  g/mol 

32  g/mol 

cd 

0.967  ±0.048 

1.032  ±0.052 

1.024  ±0.051 

Po 

46.70  ±  0.12  psia 
3.220E5  ±  800  N/m2 

124.38  ±0.31  psia 
8.576E5±  2100  N/m2 

68.55  ±  0.17  psia 
4.726E5  ±  1200  N/m2 

To 

294  ±  3  K 

292  ±  3  K 

295  ±  3  K 

8N/  a  r 
/8Cd 

5.53  mmol/s 

25.55  mmoEs 

5.33  mmol/s 

8N/  a p 
/8P0  0 

0.28  mmol/s 

1.28  mmol/s 

0.27  mmol/s 

dN/  at 
/8T0  0 

-0.55  mmol/s 

-2.56  mmol/s 

-0.53  mmol/s 

N 

110.7  ±  5.6  mmol/s 

511.1  ±  25.7  mmol/s 

106.6  ±  5.4  mmol/s 

AN/ 

Zn 

5% 

5% 

5% 

A.4  Iodine  Flow  Rate  and  its  Uncertainty 

The  iodine  molar  flow  rate  is  determined  by  the  molar  flow  rate  of  the  carrier  helium  and 
the  concentration  of  iodine  in  the  carrier  flow.  The  concentration  of  iodine  is  measured  using  the 
absorbance  of  an  argon-ion  laser  beam  that  passes  through  the  flow  of  carrier  helium  and  iodine. 
This  method  is  detailed  and  patented  by  Davis.'  The  concentration  measurement  and 
uncertainty  in  the  concentration  and  iodine  flow  rate  will  be  described  below. 
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A.4.1 


Iodine  Concentration  Measurement 


The  iodine  concentration  is  measured  using  the  absorbance  of  an  argon-ion  laser  that 
passes  through  the  helium  carrier  and  iodine  flow,  using  Beer’s  Law.  A  form  of  Beer’s  Law 
is 

IT=I0e-NoL  (A.  17) 

N  is  the  number  density  of  absorbing  iodine  molecules,  a  is  the  absorbance  cross  section  of 
iodine,  L  is  the  path  length  of  the  laser  beam  through  the  flow,  and  Id  and  It  are  the  laser 
intensities  before  and  after  passing  through  the  flow.  At  488  nm,  the  absorption  cross  section 
was  measured  to  be  1.70  x  10-18  cm2  by  Davis  and  1.68  x  10‘18  cm2  by  Saiz-Lopez  et  al. 53-54 
Beer’s  law,  rearranged  to  give  the  number  density,  N,  is 

N  =  —In—  (A.  18) 

aL  IT 


The  concentration  of  iodine  in  the  flow  is  the  number  density  divided  by  Avogadro’s 
constant,  Na, 


°lna 


lo 

1T 


(A.  19) 


Figure  A.2  is  a  schematic  of  the  diagnostic  setup.  The  argon-ion  laser  beam  exits  the 
laser,  passes  through  a  100-Hz  optical  chopper,  and  a  beam  splitter  sends  50%  of  the  beam  to 
the  lo  photometer  and  the  other  50%  is  coupled  into  a  fiber  optic  cable.  The  fiber  optic  cable 
is  routed  to  the  iodine  and  helium  pipe,  where  a  cell  with  windows  is  attached.  The  laser 
beam  passes  through  the  windowed  cell  (which  has  a  path  length  of  1 .27  cm)  and  is  coupled 
into  another  fiber  optic  cable.  The  beam  is  then  directed  to  the  It  photometer.  The  signals 
from  the  photometers  are  amplified  by  lock-in  amplifiers.  A  baseline  reading  of  ItIIo  is  taken 
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with  no  iodine  present  in  the  flow.  The  ratio  is  multiplied  by  a  factor  to  result  in  a  ratio  of  1 . 
This  factor  is  used  in  the  measurements  when  iodine  is  present. 


Figure  A.2.  A  schematic  of  the  diagnostic  setup  to  measure  iodine  concentration  in  the  iodine/helium 
flow. 


A.4.2  Iodine  Concentration  Measurement  Uncertainty 

The  uncertainty  of  the  iodine  concentration  is  determined  by  the  uncertainty  of  the 
photometer  measurements  and  the  path  length  of  the  cell.  Iodine  uncertainty  is 


ay 

V  5/o 


AL 


+ 


ay 

dl-r 


Mr 


+ 


dL 


(A. 20) 


where 

ah]  i 

dIT  oLNaI0 

ah]  -i 

dl0  oLN  AIT 


(A. 21) 


(A. 22) 
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(A.23) 


ffe]  -1  lnh 

dL  oL2Na  It 

Based  on  the  measurement  standard  of  deviation,  the  uncertainty  in  the  photometer 
signal  is  ±  1%.  With  fabrication  tolerances,  the  uncertainty  in  the  cell  path  length  is  ±  0.02 
cm.  Table  A. 2  displays  two  samples  of  measured  and  calculated  values  used  in  determining 
the  iodine  concentration  and  its  uncertainty.  There  is  more  iodine  concentration  in  the 
helium  carrier  flow  with  Sample  Data  2,  which  has  significantly  less  uncertainty.  With  a 
higher  concentration  of  iodine,  the  laser  intensity  after  passing  through  the  flow,  It,  is 
reduced,  reducing  the  uncertainty  in  IT.  It  is  therefore  advantageous  to  have  a  higher 
concentration  of  iodine  in  the  flow  to  reduce  the  concentration  uncertainty. 


Table  A.2.  Sample  calculation  of  iodine  concentration  and  uncertainties. 


Sample  Data  1 

Sample  Data  2 

It 

100  ±  1 

100+1 

IT 

97.0  +  1.0 

19.0  +  0.2 

L 

1.27  ±  0.02  cm 

1.27  ±  0.02  cm 

7.69 

7.69 

-7.93 

-40.6 

dVaL 

-1.17 

-20.1 

hi 

23.4  +  11.1  mmol/m3 

1280  ±  50  mmol/m3 

47% 

4% 
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A.4.3 


Iodine  Molar  Flow  Rate  Measurement 


The  iodine  molar  flow  rate  is  determined  by  Dalton’s  Law  of  Partial  Pressures, 
definition  of  mole  fraction,  and  the  ideal  gas  law.  Per  Dalton’s  Law  of  Partial  Pressures, 
Protai  =  P,2  +  PHe  (A.  24) 

The  mole  fraction  of  iodine  is  then 


h 


(A. 25) 


The  total  pressure  is  measured  in  the  iodine  concentration  diagnostic  cell,  but  the 
partial  pressure  of  iodine  cannot  be  measured.  However,  the  ideal  gas  law  can  be  used  to 
represent  the  partial  pressure  as  a  function  of  temperature  and  iodine  concentration.  The 
temperature  is  measured  with  a  thermocouple  and  the  iodine  concentration  is  calculated  from 
the  argon-ion  laser  beam  absorbance, 

Ph  =RTn/V  =  RT[l2]  (A. 26) 


Therefore,  the  iodine  mole  fraction  is 


x,  = 

17 


RT[lA 


Total 


(A. 27) 


Now  that  the  iodine  mole  fraction  can  be  calculated,  the  iodine  molar  flow  rate  must 
be  detennined.  The  iodine  mole  fraction  is  also  a  function  of  the  helium  and  iodine  flow 
rates,  where  A  is  a  molar  flow  rate 


N, 

17 


* /,  +  *«. 


(A. 28) 


Equation  A.28  can  be  solved  for  iodine  molar  flow  rate  as  a  function  of  the  helium  molar 
flow  rate  and  iodine  mole  fraction 
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N,  =N 

17 


He 


1  -  x. 


(A. 29) 


l2  J 


Inserting  Eqn.  A. 27  into  Eqn.  A.29  for  mole  fraction  gives 


Nr  =  NHe 


r  pTotaI  ,vl 

[I2]RT 


(A. 30) 


Equation  A. 30  gives  the  iodine  molar  flow  rate  as  a  function  of  the  measured  pressure, 
temperature,  and  iodine  concentration  in  the  diagnostic  cell  and  the  helium  molar  flow  rate. 


A.4.4  Iodine  Molar  Flow  Rate  Measurement  Uncertainty 

The  uncertainty  of  the  iodine  molar  flow  rate  is  determined  by  the  uncertainties  of  the 
iodine  concentration,  measured  pressure  and  temperature  in  the  iodine  diagnostic  cell,  and 
helium  flow  rate.  Iodine  molar  flow  rate  uncertainty  is 


A Nj  = 


dN,  r  , 


dN, 


dR 


-A R 


Total 


Total 


+ 


dN, 

1 2 

8T 


-AT 


dN, 


dN 


-AN 


He 


He 


(A. 31) 


where 


8Nh  _  »v  PT„„, 

SIR  "■ 


;[i2]2rt 


(A.32) 


dN, 

_ h_ 

dR 


=  -N 


i/[i2]rt 


He 


Total 


(A. 33) 


dN, 

12 

dT 


=  N 


Total 


[r_]rt2 


He 


(A. 34) 
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1 


(A. 35) 


dN, 

l2  _ 

BN„,  ~  PMllh¥T~  1 

Table  A. 3  displays  the  measured  and  calculated  values,  along  with  their  uncertainties, 
used  in  calculating  the  iodine  molar  flow  rate  and  its  uncertainty.  Two  sample  calculations 
were  done,  both  of  which  correspond  to  the  sample  calculations  from  Table  A. 2.  The 
calculation  from  Sample  Data  2  has  a  significantly  smaller  uncertainty  percentage  than  that 
from  Sample  Data  1 .  This  is  because  the  percentage  of  uncertainty  in  the  iodine 
concentration  is  significantly  smaller  in  Sample  Data  2. 


Table  A.3.  Two  sample  calculations  of  iodine  molar  flow  rate  and  their  uncertainties. 


Sample  Data  1 

Sample  Data  2 

hi 

23.4  ±  11.1  mmol/m3 

1280  ±  50  mmol/m3 

P  Total 

7.82  ±  0.02  psia 

53900  ±  140  N/m2 

13.90  ±  0.03  psia 
95800  ±  240  N/m2 

T 

402  ±4 

429  ±4 

NHe 

1 10.7  ±  5.6  mmol/s 

79.7  ±  4.0  mmol/s 

/4[/2  ]Ah= i 

7.60  x  10"3  mmol/s 

1.50  x  10"1  mmol/s 

A/  a P 

/ ^  Total 
/  Total 

4.03  x  10"4  mmol/s 

1.04  x  10"2  mmol/s 

8Nh/mKF 

1.61  x  10"3  mmol/s 

4.18  x  10"2  mmol/s 

ax/4> 

8.09  x  10"3  mmol/s 

2.00  x  10"1  mmol/s 

Nj 

l2 

0.16  ±  0.08  mmol/s 

3.98  ±  0.25  mmol/s 

AN,  / 

A, 

50% 

6% 
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A.5  Injector  Molar  Flow  Rate 


To  accurately  deliver  a  relatively  small  quantity  of  iodine  and  carrier  helium,  it  is  best  to 
flow  a  large  quantity  (at  the  desired  ratio  of  iodine  to  helium)  and  fractionate  the  flow  so  that  the 
desired  fraction  of  this  flow  is  delivered  to  the  injectors  whereas  the  rest  is  passed  through  to  the 
vacuum.  The  iodine  and  carrier  helium  flow  is  delivered  at  an  approximate  rate  of  125  mmol/s 
of  helium  and  1 .9  mmol/s  of  iodine.  For  a  single  injector,  for  example,  only  1%  of  the  flow 
should  be  delivered  to  the  injector,  while  99%  of  the  flow  is  sent  to  the  vacuum.  As  seen  in  Eqn. 
A.  10,  the  flow  rate  is  dependent  on  the  area  of  an  orifice.  In  the  situation  where  flow  needs  to  be 
fractionated,  two  orifice  areas,  one  area  from  the  injector(s)  and  another  from  the  bypass  orifice, 
make  up  the  total  area  used  in  Eqn.  A.  10.  The  molar  flow  rate  based  on  the  sum  of  a  single 
bypass  orifice,  Abp,  and  the  injector  throat  area,  Ainj,  is 


N  =  CdP0(Ahp+AinJ). 


y  f  y  + 1  'j  /2^~y) 


HIM  T0  2 


(A. 36) 


For  simplicity,  this  equation  is  simplified  with  a  constant,  B  ’ 


»  =  ^B'(Abp+AJ 


(A. 37) 


where 


y  f  y  + 1  'j  /^(wO 


9tM\  2 


(A. 38) 


The  ratio  of  injected  molar  flow  to  the  total  molar  flow  is 


C  P 
AT  B' A 


N  CdPn 


BiAhP+AJ  bp 


+  A:, 


(A. 39) 
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With  area  directly  proportional  to  the  square  of  the  orifice  (or  injector)  throat  diameter,  the 
injector  molar  flow  rate  is  written  as 


N 


d-, 


dlp 


+  d  ■, 


(A. 40) 


where  dmj  and  dbP  are  the  injectors’  effective  throat  diameter  and  the  bypass  throat  diameter, 
respectively. 


A.6  Injector  Molar  Flow  Rate  Uncertainty 

The  uncertainty  in  the  injector  and  bypass  diameters  and 


and  iodine  molar  flow  rate  dictate  the  uncertainty  in  the  injector 


for  the  injector  molar  flow  rate  is 


AN.  .  = 

inj 


rdN. . 

— —A  d. 

yddinj 


rdN.  . 

— —Ad 

yddbp 


bP 


+ 


dN„ 


inj 


dN 


AN 


where 


dNinj  _ 

2dbpdinj 

dd  ini  ( 

' d l  +J2.)2 

v  bp  in]  ) 

dN.  . 

2d 2  d . 

in] 

inj  bp 

ddbp  “  j 

[dl+dl] 

dN.  ■ 

d2. 

in] 

inj 

dN  d2p  +  d;nj 


bN  =  +{&Nj 


the  uncertainty  in  the  helium 
flow  rate,  Ninj .  The  uncertainty 

(A.41) 

(A.42) 

(A.43) 

(A. 44) 

(A.45) 
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Sample  Data  2  did  not  fractionate  the  flow,  but  Sample  Data  1  did.  Table  A.4  displays  values 
used  in  calculating  the  injector  flow  rate.  The  uncertainty  of  5%  is  relatively  small. 

Table  A.4.  Sample  calculation  of  an  injector  molar  flow  rate  and  its  uncertainty. 


Sample  Data  1 

dinj 

0.020  ±0.001in 

dbp 

0.160  ±0.001in 

N 

1 10.9  ±  5.6  mmol/s 

'S’ 

<1 

'S’ 

1.51  x  10"3  mmol/s 

£ 

^3 

<1 

.§• 

\ 

\(0 

=A. 

ro 

1.89  x  10"4  mmol/s 

dl*iny(-AN 
/  dN 

8.57  x  10"2  mmol/s 

Nmj 

1.71  ±  0.09  mmol/s 

AAA/ 

/Ninj 

5% 

A.7  Injector  Mole  Fraction  Uncertainty 

Now  that  the  total  flow  rate  through  the  injector  is  known,  the  mole  fractions  of  iodine 
and  helium,  with  the  corresponding  uncertainties  must  be  detennined.  The  mole  fractions  are 
calculated  with  the  molar  flow  rates.  The  mole  fractions  for  helium  and  iodine,  respectively,  are 


N 


XHe  = 


He 


NHe  +  N, 


(A. 46) 


xi 


NHe  +  N, 


(A. 47) 
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The  uncertainties  for  the  mole  fractions  are 


Ax,  = 
l2 


where 


dXHe  _ 

dxh 

SNHe 

(A  +  A) 

2“  snH( 

OX  He  _ 

-A« 

dx, 

_  l2 

ON, 

12 

(A,  +  A  )2 

dN, 

12 

(A. 49) 

(A. 50) 

(A. 51) 
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Table  A.5.  Sample  calculation  of  the  injector  mole  fractions  and  their  uncertainties. 


Sample  Data  1 

Nj 

l2 

0.16  ±  0.08  mmol/s 

NHe 

1 10.7  ±  5.6  mmol/s 

dXHe/ . 

/dNHe  He 

7.29  x  10-5 

/dN,  h 

-6.89  x  10"5 

AN 

/snh,  “ 

-7.29  x  10"5 

&'/4  AW, 

/  <9tV,  - 

/  12 

6.89  x  10-5 

XHe 

0.99858  ±  0.00069 

XI 2 

0.00076  ±0.00069 
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Appendix  B:  Flow  Straightener 


B.l  Introduction 

The  primary  flow  delivered  to  the  nozzle  is  a  mixture  of  helium  and  oxygen.  To 
adequately  mix  the  two  gases,  an  in-line  mixer  is  used  to  induce  swirl  in  the  flow.  Once  the  two 
gases  are  mixed,  the  turbulence  generated  by  the  in-line  mixer  is  no  longer  desired.  Mesh  wire 
screens  and  honeycomb  can  be  placed  in  the  flow  to  suppress  turbulence.  However,  care  must  be 
taken  in  choosing  the  mesh  size  of  the  screens,  honeycomb  size,  honeycomb  length,  and  relative 
placement  between  screens  and  honeycomb  because  more  turbulence  can  be  produced  than 
suppressed  if  these  design  elements  are  not  chosen  with  care. 


B.2  Flow  Straightener  Design 

The  mixture  of  gases  is  delivered  in  a  pipe  that  has  a  diameter  of  1 .0  in  to  a  plenum  that 
is  rectangular  in  shape,  with  a  width  of  1 .669  in  and  a  height  of  2.204  in.  A  transition  stage  is 
used  to  expand  the  circular  pipe  to  a  rectangular  cross  section,  matching  the  dimensions  of  the 
plenum.  The  elements  used  to  reduce  the  flow  turbulence  are  placed  between  this  transition 
stage  and  the  entrance  of  the  plenum. 

Farell  and  Youssef  recommend  placing  a  combination  of  honeycomb  and  mesh  screens  to 
reduce  both  lateral  and  axial  turbulence.5^  The  honeycomb  reduces  swirl  and  large-scale  eddy 
motions  while  the  mesh  screens  reduce  the  axial  variations  in  velocity.  In  their  study,  the  best 
combination  for  reducing  turbulence  was  a  coarse  screen,  honeycomb,  fine  mesh  screen,  and 
even  finer  mesh  screen.  In  their  study,  the  turbulence  intensities  (ratio  of  velocity  fluctuations  to 
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average  velocity)  were  reduced  from  12  to  2  percent.  They  tried  other  combinations  of 
honeycombs  and  screens,  but  this  configuration  had  the  most  significant  effect  on  decreasing 
turbulence. 

Loehrke  and  Nagib  studied  the  effect  honeycomb  size  had  on  suppressing  turbulence.56 
While  a  relatively  long  honeycomb  reduces  internal  turbulence  more  than  shorter  honeycomb, 
the  long  honeycomb  generates  turbulence  through  the  shear  layer  instability,  such  that  the  net 
reduction  of  turbulence  by  the  long  honeycomb  is  less  than  with  a  shorter  honeycomb.  They 
recommend  the  honeycomb  length  to  be  less  than  10  cell  sizes  long  (though  their  best  results  had 
a  length  to  cell  diameter  ratio  of  6).  Their  results  include  the  turbulence  decay  length  for 
honeycomb  of  various  sizes,  where  the  honeycomb  with  the  length  to  cell  ratio  of  6  had  a  decay 
length  of  about  6  cell  lengths.  They  also  recommend  using  a  mesh  screen  downstream  of  the 
honeycomb,  placed  less  than  5  cell  sizes  downstream.  Scheiman  and  Brooks  studied  honeycomb 
and  mesh  screen  combinations  for  reducing  flow  turbulence  and  conclude  that  the  honeycomb 
that  is  used  to  reduce  lateral  turbulence  should  be  followed  by  a  screen  to  reduce  the  axial 
turbulence  created  by  the  honeycomb.57  With  their  published  data,  it  is  evident  that  a  finer  mesh 
screen  has  greater  turbulence  reduction.  Two  mesh  screens  equal  in  size,  or  in  cascade  starting 
with  the  largest  mesh  first,  can  be  used  to  suppress  the  axial  turbulence.  Laws  and  Livesey 
recommend  the  second  screen  to  be  twice  as  fine.  They  also  recommend  placing  the  screens 
500  mesh  lengths  apart  to  allow  the  turbulence  from  the  first  screen  to  decay  before  approaching 
the  second  screen. 

With  the  recommendations  from  Fare  11  and  Youssef,  Loehrke  and  Nagib,  Scheiman  and 
Brooks,  and  Laws  and  Livesey,55'58  the  flow  straightening  system  was  designed  with  a 
honeycomb  (HC)  positioned  at  the  exit  of  the  transition  stage  to  take  out  the  initial  swirl  of  the 
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in-line  mixer,  a  spacer  longer  than  the  decay  rate  of  the  honeycomb  (SI),  another  honeycomb  of 
the  same  size  to  suppress  any  remaining  swirl,  a  fine  mesh  screen  (FS1),  a  spacer  of  about  500 
mesh  lengths  (S2),  and  an  even  finer  mesh  screen  (FS2).  The  second  mesh  screen  is  positioned 
at  the  entrance  of  the  plenum.  The  plenum  is  longer  than  the  decay  rate  of  the  second  mesh 
screen.  The  properties  for  the  flow  straightening  elements  are  displayed  in  Table  B.l.  Figures 
B.l  and  B.2  display  a  schematic  and  a  picture  of  the  flow-straightening  elements. 

The  effective  diameter  of  the  plenum  is  1 .90  in,  as  determined  by 

-  p  ’ 

where  A  is  the  area  and  P  is  the  perimeter  of  the  plenum.  The  honeycomb  cell  diameter  is  6.6% 
of  the  effective  plenum  diameter,  which  is  similar  to  the  5.7%  honeycomb  cell  to  pipe  diameter 
used  by  Farrell  and  Youssef.  The  cell  length  to  diameter  is  6.4,  which  is  recommended  by 
Loehrke  and  Nagib.  The  spacer,  SI,  is  6.25  honeycomb  cell  diameters  long,  which  is  the 
measured  decay  length  in  Loehrke  and  Nagib ’s  work.  The  second  fine  mesh  screen  (FS2)  is 
about  twice  as  fine  as  the  mesh  in  FS1,  as  recommended  by  Laws  and  Livesey.  In  choosing  the 
sizes  of  the  honeycomb  and  mesh  screens,  it  was  important  to  choose  elements  that  would  not 
decrease  the  plenum  area  too  much  to  have  a  significant  pressure  drop.  The  finest  mesh  screen 
(FS2)  restricts  the  area  the  most  with  an  open  area  of  36%  of  the  plenum  area,  increasing  the 
Mach  number  from  0. 12  to  0.33,  which  is  still  reasonable  for  the  upstream  portion  of  the  plenum. 
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Table  B.l.  Flow  straightening  element  properties. 


Length  (in) 

Cell/Mesh  Size  (in) 

Foil/Wire  Thickness  (in) 

HC 

0.8 

0.125 

0.004 

SI 

5.0 

N/A 

N/A 

FS1 

N/A 

0.0041 

0.0026 

S2 

2.0 

N/A 

N/A 

FS2 

N/A 

0.0024 

0.0016 

Figure  B.l.  A  schematic  of  the  flow  straightening  elements. 


Figure  B.2.  A  picture  of  the  flow  straightening  elements. 
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Appendix  C:  Nozzle  Design 


C.l  Introduction 

The  nozzle  profile  was  designed  to  have  a  contour  on  the  top  and  bottom  of  the  nozzle 
but  flat  sidewalls.  The  nozzle  profile  was  calculated  and  then  corrected  to  account  for  the 
boundary  layer  that  would  build  up  on  it.  In  this  appendix,  the  design  of  the  nozzle  contour  and 
boundary  layer  correction,  from  the  nozzle  throat  to  nozzle  exit  plane,  is  discussed. 


C.2  Nozzle  Design 

The  contour  of  the  nozzle  was  designed  using  a  computer  program  developed  at  the 

TO 

University  of  Illinois  at  Urbana-Champaign.  The  computer  program  is  written  in  FORTRAN 
and  implements  a  method  of  characteristics  design  analysis  for  a  supersonic  nozzle  that  has  a 
continuous  slope  at  the  throat.  Input  conditions  include  the  ratio  of  specific  heats  for  the  gas 
used  in  the  design  nozzle,  desired  Mach  number,  throat  height,  radius  of  curvature  at  the  throat, 
and  wall  shape  (hyperbolic  or  parabolic).  The  program  calculates  a  nozzle  contour  that  results  in 
a  uniform  exit  flow  aligned  with  the  nozzle  axis.  For  the  nozzle  used  in  this  research,  the  input 
conditions  are  displayed  in  Table  C.l  and  the  nozzle  contour  is  plotted  in  Figure  C.l. 
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Nozzle  Profile  (mm) 


Table  C.l.  Nozzle  design  code  inputs. 


Code  Input 

r 

1.587 

M 

2.2 

Throat  Height 

10mm 

Throat  Radius 
of  Curvature 

50mm 

Initial  Wall 
Shape 

Hyperbolic 

Distance  from  Nozzle  Throat  (mm) 
Figure  C.l.  Nozzle  contour  from  method  of  characteristics  code. 


C.3  Boundary  Layer  Correction 

To  maintain  the  design  Mach  number,  the  nozzle  contour  was  adjusted  to  account  for  a 
boundary  layer  that  would  build  up  on  the  sidewalls  and  contour.  The  boundary  layer  build-up 
was  calculated  by  estimating  the  nozzle  sidewalls  and  contour  as  flat  plates.  Using  Blasius’  flat 
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plate  theory,  the  boundary  layer  displacement  thickness,  S,  is  calculated  at  each  point  along  the 
contour  and  sidewalls  and  is  a  function  of  the  distance  from  the  throat,  x,  freestream  velocity  at  x, 
Uoo,  and  kinematic  viscosity  at  x,  k59  It  can  also  be  rewritten  as  a  function  of,  Rex,  the  Reynolds 
number  based  on  jc,  as 

<5  =  1.7208  I— =  1.7208-4=  (C.l) 

w.  ^ 

The  displacement  thickness  is  a  good  measure  for  expanding  the  contour  because  it  is  the 
displacement  thickness  of  the  external  potential  field  of  the  flow  caused  by  the  reduction  of 
velocity  in  the  boundary  layer. 

The  freestream  velocity  is  a  function  of  the  nozzle  Mach  number,  M,  and  the  speed  of 
sound,  a.  The  speed  of  sound  is  a  function  of  the  ratio  of  specific  heats,  the  gas  constant,  R,  and 
temperature 

Uao=Ma  =  MjyRT  (C.2) 


For  the  20%  oxygen  and  80%  helium  mixture,  y  is  1.587  and  R  is  866  N-m/kg-K.  The  Mach 
number  and  temperature  are  functions  of  y  and  the  area  ratio  of  the  nozzle  at  x  and  at  the  throat  (x 
=  0),  A/A*.  Isentropic  relationships  for  A/A*  and  T/To  as  a  function  of  M,  where  To  is  the  plenum 
(stagnation)  temperature,60  are 


A 

z 


M 


7  + 1  v 


J 


T 

T 

1o 


1 


r 

1  + 


y- 1 


Mz 


(C.3) 

(C.4) 
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Mach  Number 


For  this  nozzle  contour,  the  Mach  number  as  a  function  of  x  is  plotted  in  Figure  C.2  and  the 
temperature  (with  a  plenum  temperature  of  300K)  is  plotted  in  Figure  C.3.  The  velocity  (as 
calculated  in  Eqn.  C.2)  is  plotted  in  Figure  C.4. 


l  _ l _ l _ l _ l _ l _ l _ l _ l _ 

0  5  10  15  20  25  30  35  40 

Distance  from  Nozzle  Throat  (mm) 


Figure  C.2.  Mach  number  versus  distance  from  nozzle  throat. 


45 


147 


Velocity  (m/s)  Temperature  (K) 


240 

220 

200 

180 

160 

140 

120 

0  5  10  15  20  25  30  35  40  45 

Distance  from  Nozzle  Throat  (mm) 

Figure  C.3.  Temperature  versus  distance  from  nozzle  throat,  with  plenum  temperature  of  300  K. 
950 

900 

850 
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Distance  from  Nozzle  Throat  (mm) 

Figure  C.4.  Freestream  velocity  versus  distance  from  nozzle  throat. 
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The  final  variable  needed  to  calculate  the  displacement  thickness  is  the  kinematic 


viscosity,  v 


(C.5) 


which  is  viscosity,  //,  divided  by  the  gas  density,  p.  The  density  is  calculated  as  a  function  of  x 
using  the  Ideal  Gas  Law 


P 

RT 


(C.6) 


where  pressure,  P  was  determined  using  the  isentropic  relationship  between  the  pressure, 
stagnation  pressure,  and  Mach  number. 


(C.7) 


The  plenum  pressure  is  65  Torr  (8666  Pa).  Plots  of  the  pressure  and  density  are  shown  in 
Figures.  C.5  and  C.6,  respectively. 
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Density  (kg/m  )  Pressure  (Pa) 


4500 

4000 

3500 

3000 

2500 

2000 

1500 

1000 

500 

0  5  10  15  20  25  30  35  40  45 

Distance  from  Nozzle  Throat  (mm) 

Figure  C.5.  Freestream  pressure  versus  distance  from  nozzle  throat. 
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Figure  C.6.  Freestream  density  versus  distance  from  nozzle  throat. 
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The  viscosity  is  calculated  for  the  mixture  of  oxygen  and  helium.  Using  Sutherland’s 
Law,  the  viscosity  for  each  gas  is  determined  as  a  function  of  temperature  and  the  Sutherland 
constant,  S 61 


JL 

Mo 


3/2 


V^o  J 


T0+S 
T  +  S 


(C.8) 


For  oxygen  and  helium,  the  Sutherland-law  viscosity  parameters  are  displayed  in  Table  C.2.  The 
viscosities  for  each  gas  are  plotted  in  Figure.  C.7. 


Table  C.2.  Sutherland-law  viscosity  parameters.62 


Helium 

Oxygen 

To  (K) 

273 

273 

jUo  (N-s/m2) 

2.04E-5 

1.92E-5 

S(  K) 

36 

139 

151 


X  10 

2  - 1 - 1 - 1 - 1 - 1 - 1 - 1 - r 
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Figure  C.7.  Viscosity  for  the  oxygen  and  helium  gases  versus  distance  from  nozzle  throat. 


Wilke’s  formula  for  calculating  the  viscosity  for  a  mixture  of  gases  was  used,  where  n  is  the 
number  of  gases  present  and  i,j  of  1  represents  oxygen  and  i,j  of  2  represents  helium61 


sh  xiM, 

,=1 2>M/ 

j= i 


(C.9) 


where 


(8  +  SMi/MJf2 


(C.10) 


For  the  situation  when  i  =j,  Eqn.  C.  10  reduces  to  </>=  \.  The  values  for  (/)n  and  foi  are  plotted  in 
Figure  C.8  and  the  viscosity  for  the  oxygen  and  helium  mixture  is  plotted  in  Figure  C.9.  The 
kinematic  viscosity  is  plotted  in  Figure  C.10. 
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Viscosity  (N-s/m“) 
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Figure  C.8.  Values  for  ( i>i2  and  </ h?  versus  distance  from  nozzle  throat. 
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Figure  C.9.  Viscosity  for  the  helium/oxygen  mixture  versus  distance  from  nozzle  throat. 
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Kinematic  Viscosity  (nT /s) 
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Figure  C.10.  Kinematic  viscosity  for  the  helium/oxygen  mixture  versus  distance  from  nozzle  throat. 

With  the  freestream  velocity  and  kinematic  viscosity,  the  Reynolds  number  based  on  the 
distance  from  the  nozzle  throat,  x,  can  be  calculated,  as  shown  in  Figure  C.  1 1 . 
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Figure  C.ll.  Re.t  versus  distance  from  nozzle  throat. 

From  Eqn.  C.  1,  the  displacement  thickness  can  now  be  determined,  which  is  plotted  in  Figure 
C.  12.  Figure  C.  13  displays  the  original  nozzle  profile  (as  plotted  in  Figure  C.  1)  and  the  nozzle 
profile  that  is  corrected  by  the  boundary  layer  displacement  thickness. 
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Figure  C.12.  Displacement  thickness  versus  distance  from  nozzle  throat. 
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Figure  C.13.  Original  nozzle  profile  and  nozzle  profile  corrected  by  the  displacement  thickness  of  the 
boundary  layer. 
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